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HOT-WIRE CALIBRATION IN SUBSONIC/TRANSONIC FLOW REGIMES 


by 

K.A. Nagabushana and R.L. Ash 
Old Dominion University 
Department of Aerospace Engineering 
Norfolk, VA 23529 


Abstract ^ 

A different approach for calibrating hot-wires, which simplifies the calibration procedure 
and reduces the tunnel run-time by an order of magnitude was sought. In general, it is accepted 
that the directly measurable quantities in any flow are velocity, density, and total temperature. 
Very few facilities have the capability of varying the total temperature over an adequate range. 
However, if overheat temperature parameter, a w , is used to calibrate the hot-wire then the 

directly measurable quantity, voltage, will be a function of the flow variables and the overheat 
parameter i.e., E- f{u,p,a w ,T w ) where a w will contain the needed total temperature 

information. In this report, various methods of evaluating sensitivities with different dependent 
and independent variables to calibrate a 3-wire hot-wire probe using a constant temperature 
anemometer (CTA) in subsonic/transonic flow regimes is presented. The advantage of using a w 
as the independent variable instead of total temperature, T 0 , or overheat temperature parameter, 
t , is that while running a calibration test it is not necessary to know the recovery factor, the 

coefficients in a wire resistance to temperature relationship for a given probe. It was deduced 
that the method employing the relationship E = f{u,p,aj) should result in the most accurate 

calibration of hot wire probes. Any other method would require additional measurements. Also 
this method will allow calibration and determination of accurate temperature fluctuation 
information even in atmospheric wind tunnels where there is no ability to obtain any temperature 
sensitivity information at present. This technique greatly simplifies the calibration process for 
hot-wires, provides the required calibration information needed in obtaining temperature 

fluctuations, and reduces both the tunnel run-time and the test matrix required to calibrate hot- 

/ 

wires. Some-of-the-results using the above techniques are presented in an appendix. 


Nomenclature 


R 

a — — , Overheat temperature parameter 

* 

aaw 

d w Wire diameter 

E Mean voltage across wire 

Kn Knudsen number 

k, Thermal conductivity of air at T 0 

t Wire length 

M Mach number 

m Mass of fluid flow 

Nu, Nusselt number 

P w Power to hot-wire 

Q w Heat transfer rate 

R nf Resistance of wire at reference temperature 

Re t Reynolds number based on viscosity evaluated at T a and wire diameter 
R w Resistance of wire 

S u Velocity sensitivity 

S p Density sensitivity 

S To Total temperature sensitivity 

Recovery temperature of wire 
T 0 Total temperature 

T„f Reference temperature 

T w Wire temperature 

u Velocity 


a&P Coefficients in wire resistance to temperature relationship 



V 

p 

Y 
0 

T 



— , Recovery temperature ratio 


Density 

Ratio of specific heats 
T w 

— , Temperature parameter 
T -T 

w ~ ~ ■ Temperature parameter 

o 


Introduction 


The heat transfer from heated wires in fluid flows has been investigated for many years 
for various reasons. One of the main reasons in fluid mechanics is the application of the results 
to hot wire anemometry. The early work in this field was performed by Boussinesq (1905) 1 and 
King (1914) 2 . Even though " King's /aw" was derived over 75 years ago, it is still being used, in 
various forms, to predict the heat transfer from heated wires in subsonic incompressible flows. 
From King's law K can be shown that the heat transfer parameter, Nusselt number, is a function 
only of Reynolds and Prandtl numbers. 

In the 1950's, it was shown that the Nusselt number can be a function of Mach number 
for Mach numbers as low as 0.10, indicating that King's law does not apply 3 ' 4 . The reason for 
this disagreement was determined by Spangenberg 5 to be due to gas rarefaction effects. 
Spangenberg found that Nu t = f(Re t ,M,Kn), where the apparent Mach number or 

compressibility effect was a slip flow or Knudsen number effect. Since the Knudsen number is 
related to the Reynolds and Mach numbers, only two of the three quantities in the relationship for 
the Nusselt number are required. Experimental results obtained by King and others indicate that 

the Nusselt number is also a function of a temperature parameter. The functional relationships 
expressed for Nusselt number are Nu t = or Nu t = f(M,Kn, r) 7 A From the 


3 



above functional relationship, it has been shown that the mean voltage measured across the 
heated wire is E- f(u,p,T a ,T w ) 6 . 


Theoretical Consideration 

To obtain the sensitivity relationships for use in hot wire anemometry, the following heat 
transfer equation is usually used to describe the heat loss from a heated wire mounted normal to 
the flow: 


Q w = 27r£k t (T w -rjT o )Nu t 


(D 


For steady state conditions the electrical power supplied to the wire is equal to the heat 
loss from the wire (P w = Q w ). Differentiating this equation in terms of its logarithmic derivatives 

yields: 


d log Q w = d\ogNu t ~—d\ogT o + —d\og 7^ - — d\og tj+ d\og k t (2) 

T T 7 


where d\ogT w =0 for a constant temperature anemometer (CTA). 
For sensitivities depending on Nu ( = t): 


Equations were re-derived for the sensitivities of heated wires operated in compressible 
flow and powered by a constant temperature anemometer using the heat transfer relationship for 
the variables Nu t = f(M,Kn, r). The sensitivity equations derived by Baldwin? and 

Sandbom 8 are incomplete for several reasons. In these references the velocity sensitivity 
includes a cTTjdu term under the constraint that T a = constant. There are other terms missing 
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in the equations due to an incomplete partial differentiation process. Following the 
assumptions?- 8 that Nu t = f(M,Kn,r), the change in Nusselt number is: 


dlog Nu t 


dlogM+ ^k^ t n ogKn+ ^3^ dlog T 

dlog M dlog Kn c\ og x 


(3) 


The Mach number is a function of u and T c . Thus: 

dlog M = ^L d X og u + ^rd\og T a (4) 

dlogU dlog T a 

JlogM = fl+^M 2 ^logu-|fl+^M 2 \/logr o (5) 


The Knudsen number is a function of only density for a given wire diameter, and can be 
expressed as?: 


Kn = 


1. 587x1 0~ 8 

pd w 


or d\o%Kn = -d\ogp 


( 6 ) 


The recovery temperature ratio, rj , is assumed to be a function of M and Kn. Then: 


dlog rj = d\og M + -——^—^—dlog Kn 

5 dlog M 6 dlog Kn 6 


(?) 


The thermal conductivity for air is evaluated at T , giving: 


di 0g k, 

dlogT a 


( 8 ) 
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The power to the heated wire can be written as: 

d\og P w = d\og E w - d\og (9) 

For a constant temperature anemometer: 

r/log/^ = 0 

T -tjT 

The quantity x= — — — therefore, the change in x is: 

o 

d\og t= ~—d\og T ~—d\og T] 

X X 

Equations (5-1 1) can be substituted into (2) to give: 

d\ogE w = S u d\og u+Spd\ogp+S T ' d\og T a (12) 


( 10 ) 


( 11 ) 


where 



2 



M 1 


dlog Nu, 
c\ogM 


rj dlog 1 1 ' d\agNu t + ^ 
x dlog M \ diog x ) 


(13) 


s 

p 


l_ 

2 


dlogNu t ' r] dlogrj r dlogNu, 
dlog Kn x dlog Kn ^ dlog x 



(14) 



rj | dlogk, 

2[ 2 V rj dlogx x dlogT a 


(15) 
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The sensitivity equations as derived by Baldwin in reference 7 using functional 
relationship Nu, = r) are as follows: 


<n_ 

i 


1 




c\o%Nu, dp 
dlog Kn p 


r d\o% Nu, 
< cAogM 

r cfto%k t 
^Og T a ~ 


)m 2 

X 34) X 


tj x+ij (Nu, ^ dT 0 
x Nu t 34 J T a 


dU 

U 


(16) 


These equations (13-15) for velocity, density and total temperature sensitivities are 
different from equation (16) for several reasons. In equation (16) the velocity was differentiated 
with respect to T a under the assumption that T a = constant. Also, the differentiation processes 

in was incomplete. 


The re-derived equations were used to obtain the velocity, density and total temperature 
sensitivities for hot-wires and are compared with the sensitivities computed using the functional 
relationship E = f(u,p,T 0 ,T w ). The constraints required by calculus and the functional form of 

the equation dictated by the data were met by both methods for computing sensitivities. These 
results are presented in the enclosed appendix. 

For sensitivities depending on Nu, = f(M,Kn,a w ): 


Following the assumptions that Nu, = f(M,Kn,a w ), the change in Nusselt number is: 


dlogNu, = a ° gN ^ d\ogM+ d\og Kn + - ^° g Nu ‘ - d\oga w (17) 

cAog M d log Kn a log a w 


R — R 

The quantity a w = — — and = f(M,Kn,T a ). Therefore, the change in a w is: 

R a^ 
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ftogR 


j gdw 1 


<floga w = 


9 


cAogM 


\ + zzl^ 


</logu-i-i#log3; 

l £ 




a ° S -^ dlogp+ ^° gR - ^ d\ogT 0 


d\og Kn 


&0gT o 


(18) 


Substituting equations (5-10), (17) and (18) in equation (2) leads to equation (12) where 


s --\{' +z t m2 )_ 


c\ogNu t T] diogrj dlogR^ 6 c\ogNu t 

c\ogM x cAogM diogM a w d\ogM J 


(19) 


s _ 1 f ftogNu t | TJ dlog T] | 0 dlog Nu, dlogR^ 
P 2 [ <51og A>i t S\ogKn a w c\og a w PXogKn 


(20) 


= i/_2,S 0 ftogM*, 

2 } a w ftog a w x c\ogT 0 \ 


( 21 ) 


For sensitivities depending on E = f(M,Kn,a w ): 


If we use the relationship E = f(M,Kn,a w ) then the change in voltage is: 


d\ogE = 


<?log£ 

cftogM 


d\ogM + 


c\ogE 

diogKn 


d\ogKn + 


dXogE 

&oga w 


d\oga w 


( 22 ) 


Substituting equations (5) through (10) and (18) in equation (22) leads to equation (12) where 




IzIm 2 1 ^ og£ 

2 J dlog M 


9 dlogE ftogR^ 
a w c\oga w chogM 


(23) 
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s 

p 



^log E | 6 dlogE dlogR^ 
diogKn a w c\og a w c\ogKn 


( 24 ) 



0 dlogE dlogR^ } 
a w c\oga w c\ogT 0 J 


(25) 


Resistance to temperature relationship can be written as: 

R, = Wl + a(r„ - rj+/<r„- rj 2 ■ 


(26) 


Using equation (26) the quantity dlogR^j dlogT a can be obtained as: 
chogR^ T ad ^a+2fAj' adw - 7^-)j 

* ogT ° i+^-tJ+(t^-tJ 2 


For sensitivities depending on E- f(u,p,T 0 ): 


For a constant temperature anemometer (CTA), the change in voltage across the wire 
can be expressed as E = f{u,p,T 0 ). Once the wire is calibrated and functional relationships 

are defined, the sensitivities can be evaluated by taking the partial change in E with respect to 
u , p and T o in turn. The equations for sensitivities can be written as: 




U )p,T.J , 


(28) 



r dlo g£ "| 
<^Og PJu,T..T 


(29) 
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( 20 ) 


S T 


r c\ogE 
,^og T 0 


\ 

'U.p.T. 


This method of obtaining sensitivities is the direct method since the output voltage, 
which is a directly measurable quantity, is a function of velocity, density and total temperature 
which are in turn directly measurable quantities. 


For sensitivities depending on Nu t = f(M x ,Re t ,d): 


For the sake of completeness, the sensitivity equations for constant current 
anemometers (CCA) using the functional relationship of Nu, = as derived by 

Morkovin 5 , are repeated here. They are: 



f 1 | 1 »/2^ 

dlog Nu, 


V ftogT] 


4 - S 

l 2 ) 

» 

c\ogM 

R',.e 

T w eftogM 

M.O_ 

j 


1 

■■ 

<51og Nu, 


T] dlog TJ 

' 

2 

dlog Re, 

M.0 

t w <?logRe, 

M.0. 


(31) 


(32) 


S r = — 
' 2 


S, + Sj 2 m,-l)- 


c\og Nu, 


c\og 0 


M.Rt, 


(33) 


Experimental Description : 


The data were collected using the "Probe Calibration Tunner (PCT) at NASA Langley 
research center which is an open jet tunnel in which velocity, density and total temperature can 
be controlled independently. The combination of a staged pressure system, low mass flow rates 
through the facility and the availability of large vacuum systems enables this facility to run 
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continuously. It has two interchangable, subsonic and transonic nozzles, having diameters of 
1 .50 and 2.25 inches respectively, that yields a continuous flow capability over the Mach number 
range of 0.05 to 1.0. Tunnel stagnation pressure and temperature can be varied from a 
minimum of 0.20 atmospheres to a maximum of 10 atmospheres and from 500°R to 600°R, 
respectively. This corresponds to a Reynolds number range of 1 x 10^ < Ref ft <51x10^ for 

a Mach number of 1 . The 3-wire hot-wire probe utilized three different diameters of Platinum 
coated Tungsten wires (diameters of 0.0001, 0.0002 and 0.0003 inch) and was used to obtain 
most of the calibration data. The test matrix consisted of 12 velocities, 7 densities and 5 total 
temperatures for the cases where E = f{u,p, a w ,T w ). For E = f(M,Kn ,a w ) , the test matrix 

consisted of 12 Mach numbers, 7 Knudsen numbers and 8 overheats. The lower and upper 
bounds of the test matrix for the two methods are given in Tables I and II. In this study all 
Knudsen numbers were greater than 0.01, indicating that all the data were in the slip flow regime 
based on the definition of the slip flow. 


Parameter 

Minimum 

Maximum 

Total pressure, psi 

3.40306 

46.2992 

Static pressure, psi 

2.83727 

31.5238 

Velocity, ft/sec 

113.785 

950.953 

Mach number 

0.0982 

0.941 

Total temperature, °R 

500 

600 

Static temperature, °R 

417.33 

579.73 

Reynolds number/foot 

0.3 xIO 6 

12.3x10® 

Density 

0.01875 

0.14693 

Mass flow, Ib/sec 

0.012832 

1.71468 


Table I. Minimum and Maximum Test Conditions of Test Matrix for Functional 
Relationship E = f{u,p,a w ,T w ). 
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Parameter 

Minimum 

Maximum 

1 

3.675 

29.4 

Static pressure, psi 

3.44859 

26.2216 

Mach number 

0.1 

0.95 

Total temperature, °R 

540 

540 

Static temperature, °R 

457.433 

539.73 

Reynolds number/foot 

0.5 xIO 6 

13.3x10® 

Knudsen number 

0.00494 

0.10369 

Overheat ratio 

0.3 

1.0 

Density 

0.01875 

0.13333 

Mass flow, Ib/sec 

0.013103 

1.7193 


Table II. Minimum and Maximum Test Conditions of Test Matrix for Functional 
Relationship E = f(M,Kn,a w ). 


Discussion : 


To illustrate the complexity of the sensitivity equations, consider the relationship 
between voltage across a wire mounted normal to the flow and the quantities u, p and T 0 for the 
constant temperature anemometer. Utilizing a small perturbation assumption we can write: 



£ 

T 0 


(34) 
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Evaluating the partial derivatives in the above sensitivity equations requires care when 


carrying out the mean flow calibration. For example, consider equations (31-33). The evaluation 
dlogNu, 

nf ? * I krt rtktrt lnn/4 k\ i wnnilnn D rtnlu anri Mo m iimhar onH tha tntol 


dlog Re, 


must be obtained by varying P Q only and the Mach number and the total 


temperature must be held constant. Also, the evaluation of 


dlogNu, 

dlogM 


requires that the total 


pressure be changed when the Mach number is varied in order to maintain Reynolds number 

dlog 1 ) dlog 77 

constant. Similar constraints must be observed when 2 - L - and 2 - £ - are 

^og Re, M ftogM Rti 

evaluated. Care must also be taken in evaluating the partial derivatives when using other 
dependent and independent variables. This method of obtaining sensitivities is time consuming 
when conducting experiments in large wind tunnels since the Mach number and the total 
pressures must be varied in order to maintain the prescribed non-varying independent variable at 
a constant level when the variation of the remaining dependent variables being sought. 


If Nu t = f{M x ,Re t , 0) and the above described constraint is applied, the operational 

envelope for the facility must be considered, since there is a skewing of the Nuf vs. Ref curves 
for constant Mach numbers because Re, = /(«,/?) 9 . Due to this skewing, the region over 

which the partial derivatives can be evaluated will reduce rapidly. The advantage of using the 
suggested functional relationship i.e., E b = f{m,T 0 ,T w ) or E b = f{m,a w ,T w ) can be seen by 

plots of Nuf vs. Kn for constant Mach number where the data are not as skewed and a more 
complete set of derivatives can be evaluated from a given number of data points®. This efficient 
use of data with the wire voltage correlated in terms of these primitive variables cannot be 
ignored and that will reduce the amount of data needed for calibration by an order of magnitude. 
Eventually, that attribute leads to better accuracy in the calibration process. 


The effect of wire temperature is complicated due to its nonlinear variation and because 
it is not fully understood. In the literature several methods based on various heat transfer 
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relationships, using different dependent and independent variables for a given type of 
anemometer, are available to compute hot-wire sensitivities. Among these, the most extensively 
used correlations are Nu t = f(M x ,Kn, t w ) and Nu t = All these methods 

should give identical results for the sensitivity of a given hot-wire at given test conditions. Hot- 
wire calibrations have been performed as part of the PCT facility flow quality investigation, in 
various flow regimes. At higher speeds, where compressible flow effects occur, it has been 
found that King's Law is not valid 1 ' 2 - 10 . Hence, those flow regimes have been investigated. 

In the 1950's, Kovasznay 11 - 12 extended hot-wire anemometry to compressible flows 
where it was found experimentally that in supersonic flow the heated wire was sensitive only to 
mass flow and total temperature. Kovasznay developed a graphical technique to obtain these 
fluctuations, which is used primarily in supersonic and hypersonic flows. Kovasznay's 12 
compressible flow results showed that there was a significant difference between the heat 
transfer in compressible and incompressible flows. In subsonic compressible, transonic and low 
supersonic flows, effects due to compressibility influence the heat transfer from a wire. In high 
supersonic and hypersonic flows a strong shock occurs ahead of the wire and the heat transfer 

from the wire is influenced by subsonic flow downstream of the shock front. Because of this, it 
was found experimentally that in supersonic/hypersonic flows Nu, = f(Re n &) only, and the 

heat transfer from the wire was again a function of mass flow, total temperature, and wire 
temperature. For this reason some understanding of subsonic flow while studying supersonic or 
hypersonic flow, becomes essential. 

In continuum flow the mean free path of the particles is less than the diameter of the 
wire and conventional heat transfer theories are applicable. When the diameter of the wire 
approaches a few mean free paths between the particles, the flow does not behave as a 
continuum, but exhibits some effects of the finite spacing between the particles. These effects 
have been studied 13 * 14 by assuming finite velocity and temperature jump boundary conditions. 
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This gas rarefaction regime was noted as slip flow. In free molecular flow the fluid is assumed to 
be composed of individual particles and the distance between the particles is sufficiently large 
that their impact with and reflection from a body is assumed to occur without interaction between 
the particles. 

Conclusions 

The following conclusions have been made: 

• The equation obtained for a constant temperature anemometer (CTA) based on the 
assumption that E = f{u,p,T 0 ) resulted in an equation in which neither wire length nor 

the coefficients in the wire resistance vs. temperature relationship were required. This 
should result in a more accurate method for calibrating hot-wire probes in facilities where 
total temperature can be varied. 

• The equation obtained for a constatnt current anemometer (CCA) based on the 
assumption that E = f{u,p,T 0 ,T w ) results in an equation that requires the variation of 

four independent variables and is believed to be too complicated for routine use. The 
equation also requires evaluation of /? who's elimiation was desired. 

• The equation obtained for a CCA and CTA under the assumption that 
Nu t = f(M,Kn,a w ) and T a = constant resulted in an equation where a and /? occured 

in only one term in the equation for the total temperature sensitivity. Consequently, this 
calibration of wires should be more accurate than those obtained using conventional 
methods. 
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• The sensitivities computed based on the methods presented using the data obtained in the 
" Probe Calibration Tunnef are comparable. However, differences in S T ° are more 

significant than for S u and S p . 


• In general, the velocity sensitivity ranges from 0.05 to 0.40, density sensitivity ranges from 
0.10 to 0.40 and total temperature sensitivity ranges from -0.2 to -1.0. 

• The equations (13-15) for the sensitivities of a constant temperature anemometer are 

different from those in references 7 and 8. In those references, the velocity sensitivity has 

dlogT 

a term for — - under the constraints that T„ = constant. There are other terms 

dlogu 

missing in the equations due to the incomplete partial differentiation process. 
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Abstract 

Hot-wire anemometry is a well developed technique for measuring fluctuations in 
flowing fluids where its application depends on the mean flow calibration of the heated wire. The 
conventional technique usually used to calibrate hot-wire probes is based on the determination 
of Nusselt number and wire recovery temperature ratio as functions of Reynolds number and 
Mach number. This method requires the measurement of the wire length and the coefficients in 
the temperature - resistance equation for the wire material. The accurate measurement of these 
quantities are difficult to make, particularly the wire length and the coefficient of the second 
degree term in the temperature - resistance equation. The lack of accurate values for these two 
quantities are possible sources of error in the determination of the sensitivities for the heated 
wire to heat transfer from the wire. A technique is proposed that reduces the dependency of the 
sensitivities of wires to these quantities. 


Nomenclature 



c f specific heat at constant pressure 
c. specific heat at constant volume 

d w diameter of the hot-wire 

dA w elemental surface area of the wire 
e' fluctuating voltage across the sensor 
E mean voltage across the sensor 

h heat transfer coefficient 
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I current across the sensor 

k thermal conductivity of air evaluated at 

subscript temperature 
K ftogTJftogR, 

Kn Knudsen number 

l hot-wire length 

^og/f,MogT 0 
M Mach number 

Nu Nusselt number evaluated at subscript 

temperature 

P electrical power to the hot-wire 

q dynamic pressure 
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Q forced convective heat transfer 

r radius of the hot-wire 

R resistance 

Re Reynolds number 

s distance along circumference of the wire 
S m velocity sensitivity 

S p density sensitivity 

total temperature sensitivity 
T temperature 

u velocity 


x 

a 


distance along the length of wire 

(,^f 


a, & /?, coefficients in temperature - resistance 


relationship 

£ finite circuit factor, - (<5tog I /^og E) 


t) recovery temperature ratio, 7’ <jdw /7’ <> 

0 temperature parameter, R^/R^ 

p absolute viscosity of stir evaluated at 

subscript temperature 
A sweep angle 

y specific heat ratio, c r /c v 

p density 

r. temperature loading parameter, (T w - T^)/T e 
<p angle in polar co-ordinate system 

<P normalized voltage fluctuation ratio = (%) js^ 
indicates instantaneous value 


Subscrii 


adw adiabatic wall condition 

ave average condition 

/ local value 

Ml method I referred to E = f(u,p,T„T„) 

M2 method II referred to E- f(M, Kn, a^) 

o total condition 

ref reference conditions, T 0 = 0 °C 

w wire condition 

oo free stream or static condition 


the temperature - resistance relationship for the wire 
material. The operation of wires at high temperatures 
requires a knowledge of /?, values. The calibration of 
the wire at elevated temperature places a severe 
demand on the construction of the probe to withstand 
these temperatures. Because of this, values for /?, 
are often not obtained from calibration but from 
handbooks. 

The length of the wire must also be known 
and for wires having slack it is difficult to accurately 
measure their lengths. Since the wire is heated, the 
length measured when cold will differ from the value 
when the wire is hot The approximate values of /?, 
and errors in measuring the length of the wire can 
lead to error in the calibrated values of the mean flow 
sensitivities. The mean flow calibration is also 
constrained by the average values used and the 
approximations made for the heat transfer equation. 
The use of general functioned relationships for heat 
transfer from the heated wire should lead to more 
accurate mean flow calibrations. 

Two methods are proposed for calibrating the 
sensitivities of heated wires to flow variables. One 
method eliminates the need to measure the wire 
length or a, and /?,. In the second method, the 
quantities a, and /?, affect only the total temperature 
sensitivity of the wire to a minor degree. Equations, 
based on these techniques, will be presented for 
constant temperature and constant current 
anemometers (CTA and CCA). Calibration data for a 
constant temperature anemometer will be presented 
based on these two calibration techniques. 

Theoretical Consideration 
Rational for using Functional Relationships 

Consider the heat transfer from a heated wire 
with slack as shown in figure 1 . For this problem the 
general heat transfer equation due to forced 
convection can be written in integral form as: 


Introduction 


2 = \h(T,~T oiw )dA w 


( 1 ) 


Historically, much hot-wire data were 
analyzed by using equations derived by Morkovin* to 
obtain the sensitivities of the heated wire to changes 
in velocity, density, and total temperature. This is 
particularly true for wires operating in the 
compressible flow regime. The use of these equations 
requires a knowledge of the wire length, a, and /?, in 


where 

dA = dxds = rdxdtp (2) 


2 


Using equation (2) we can write equation (1) as: 


E = A«,P,T.,T.) 


( 8 ) 


Q=2rtj ) hfj.-T^dg, 


The change in voltage with respect to the 
independent variables is: 


Here h , T„ and T^, are local values. Using 
stagnation values of these quantities to normalize the 
above equation along with the stagnation values at 
zero sweep, the above equation can be rewritten as: 

n _ ... (T r > \ r *.* < r -- r ^U h (W.J a 

Q - rlh^jT. - i v A .„ (r. - T3ZZ t: (T. - TJZ t 

where A ,T„ and T waK are f(x). and 7' - . o A are 
f(M,Kn,T 0 ,\). Mach number and Reynolds number 
are f{<p). h/h^ = f(p) = f(M„Kn,) and 

Va/^.a -0 = /(Af„,An„,A) = /W. The quantities that 
are /w are due to conduction to the supports and to 
the local sweep. This is a more complex equation 
than the usual heat transfer equation used to derive 
hot-wire equations. 

The conventional heat transfer equation due 
to forced convection can be written as: 

Q = ntdJ>„{T, M - T^J (5) 

Nondimensionalizing equation (5) similar to that done 
in equation (4) gives: 

(2= ^Aa.^.-Ua^f^(6) 

"o.A-0 V. 0 

Here equation (4) and (6) represent the same quantity 
but approached in a different way. Evaluating the 
terms in equations (4) and (6) are difficult. However, 
note that both equations can be represented in terms 
of functional relationship to basic flow quantities as: 

Q= /(u.P.T.J.) (7) 


dlogE = — — d log u + ^ dlogp+ 

dog u dog p 

dog£ „ ^ dogE ,, „ 


For a constant current anemometer the 
change in wire temperature must be considered. The 
temperature of the wire should be obtained from the 
wire resistance using the following equation: 

A-=l + a 1 (7'.-rJ + A(7.-7^) 2 (10) 

Equation (10) is required to relate r. to For 

T w = /{R.) the following equation can be written: 

d\ogT w = K rflogR. (11) 

Ohm's law gives: 

dlogE. = dlogE- dlog/ (12) 

Using equation (11) and (12), equation (9) can 
be written as: 




dlogE j, . dlogE j, . dlogE .. 

— dlogu + — d\ogp+—-^—d\ogT 0 

dog u dogp dlogT c 

The final equation for the change in E becomes: 

dlogE = S„dlogn + S,,dlogp+S,.dlog7^ (14) 

where 


which can be used to evaluate the sensitivities of the 
wire to heat transfer more accurately than equation 
(4) or (5) or (6) as far as the hot-wire application is 
concerned. 


dog E/ dog u 
dogT. 




dog E/ dog p 

i-k^ i£ ( i +e ) 

dogT. 


Derivation of Benaltivit 


auationa 


The functional relationship given by Morkovin 
that governs the voltage measured across a heated 
wire was obtained from equation (5) and results in: 


dogE/dogT c 

-*■?*£( ■ ♦«> 
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For a CTA, dT„ = 0 and the final equation for 
the change in £ is same as equation (14) where: 

S # -^3LL *.*35*; S,.*2S* (16) 

dogu p dungp 5 <3og T 0 

Equation (15) is not recommended for 
calibrating a CCA since the flow variables of u , p , T c 
and T„ must be varied. This makes the calibration 
process too lengthy. The term K, in equation (11) 
also requires values of a, and ft x . However, equation 
(14) and (16) for a CTA does not require the 
knowledge of wire length or a, and ft . 

It is not often possible to vary T„ in wind 
tunnels. To overcome this problem, the following 
functional relationship is suggested: 


E = f(M,Kn,a’ w ) 


(17) 


where T e is constant Using equation (17), the change 
in the voltage across the heated wire is: 


d\og E = d log M + 

M 

dlogKn + ■^ gE d\oga'„ 
cHog Kn <9log al 

From the definition of a' w , the change in o(, is : 
d\oga' w = -^-{dlogft, - d\ogR^) 


One possible functional relationship for is: 
*ww = AM, Kn,T c ) 

From which the change in is: 


-^ R f- \-[d\ogu- ±d\ogT. 

dlog M [a), 2 


— ogR -*- d\ogp+ ^^ dlogT. 
dlogKn dlogT a 6 * 


(18) 


(19) 


( 20 ) 


( 21 ) 


Using equation (21) in (19), substituting d\oga’„ into 
equation (18), and rearranging leads to: 


dlogE = S' y dlogu + S' p d\ogp+ S^d\ogT„ + 
a. ^og a'. 


( 22 ) 


Finally, Ohm's law can be used to express the wire 
resistance in terms of the voltage across the wire as 
follows: 


dlogf?,, = (l + £)dlog£ 


(23) 


Using equation (23) in (22), the final equation for the 
change in £ is same as equation (14) where: 


5. = ASl - 


A_ 

\ <^og£ 

0 cAogE tftogR^ 

a 

[dlogM 

a'„ choga'^ cftogM 


(24) 


s, = as ' = ^i£_l (25) 

p |a(, cAogal cftogKn dlogKn) 

S z - A Si - A \— ifegg] (26) 
* * tfto%a' w \2a <5\ogM ^log7^ ) <9logA/J 

1271 

Using equation (10) the quantity dlogR,^,/ dlogT c can 
be obtained as: 


<?l°g a i + dw ~ ^./)] 

diagT. 


1 + a t ~ + A 


(28) 


The evaluation of dlogR^/dlogT^, requires 
values for a, and ft x and the latter should be obtained 
from the temperature-resistance calibration of the 
wire. The values for a, can be easily obtained, 
however, it is often necessary to use handbook values 
for ft. In any case, the values of a, and /?, occur in 
only one term in the equation for the sensitivity of the 
wire to total temperature and any error in /?, would 
be limited to this sensitivity. The possible error in 
evaluating dlogR^f c\ogT 0 by neglecting the value for 
/?, can be obtained from: 


{ T A a ' +1 ^- T ^ aiA 

glog/U. {«.^} 01 

. ^ J o I, 


( 29 ) 
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where 



and 




Kf 


- 1 


+ T nf 


(31) 


For a CTA d\ogR^ = 0, and equations (24) - (26) are 
valid with A = 1 . 


The test matrix consisted of 12 velocities, 7 
densities and 5 total temperatures for the case where 
E = f(u,p,T t ,Tj ) . For E = f{M,Kn,al) , the test matrix 
consisted of 12 Mach numbers, 7 Knudsen numbers 
and 8 overheats. The lower and upper bounds of the 
test matrix for the two methods are given in Table I 
and II. The flow region of the present test as 
compared to earlier calibration tests are shown in 
figure 5. At selected test conditions, fluctuation data 
were obtained for future analysis. In this study, all 
the Knudsen number were greater than 0.01 
indicating that all the data were in the slip flow 
regime based on the definition of the slip flow 
boundary in reference 7. 


FgdMty and Description of Present Work 

The " Probe Calibration TurtneT (PCT)^ is an 
open jet tunnel with the ability to independently 
control velocity, density and total temperature. The 
primary purpose of this facility is to economically 
calibrate probes for NASA's major test facilities at the 
Langley Research Center. Typical operational cost for 
the PCT is 1% of the cost of operating a major 
facility^. The combination of the staged pressure 
system and low mass flow requirement of the facility 
and the large volume vacuum systems enables 
economic and accurate calibration of probes. Due to 
these feature, the flow in this facility is continuous. 
Schematic diagram of the facility is presented in 
figure 2. 

The PCT has two interchangeable, subsonic 
and transonic, nozzles having diameter of 1.50 and 
2.25 inches, that gives a continuous flow capability 
over Mach number range of 0.05 to 1.0. Tunnel 
stagnation pressure and temperature can be varied 
from a minimum of 0.20 atmosphere to a maximum 
of 10 atmospheres and 500°R to 600®R, respectively. 
This corresponds to a Reynolds number range of 
1x10** < Re I ft < 51x10** for a Mach number of 1. 

A schematic diagram of the data acquisition 
system for this test is presented in figure 3. A three 
wire hot-wire probe consisting of three different 
diameters of Platinum coated Tungsten wire 
(diameters of 0.0001, 0.0002 and 0.0003 inch) was 
used to obtain the data. The 1.50 inch diameter 
nozzle was used in this test The operational envelope 
of the tunnel with this nozzle is compared with few 
other facilities in figure 4, where similar^'® hot-wire 
measurements were performed at NASA Langley. 


Parameter 


Maximum 

Total Pressure, psi 

3.40306 

46.2992 

Static Pressure, psi 

2.83727 

31.5238 

Velocity, ft/ sec 

113.785 

950.953 

Mach Number 

0.0982 

0.941 

q , psi 

0.0064 

14.328 

Total Temperature, °R 

500 

580 

Static Temperature, °R 

417.33 

579.73 

Reynolds number/ foot 

0.3 x 10 6 

12.3 x 10 b 

Heater Power, Watts 

32.5 

4545 

Density 

0.01875 

0.14693 

Mass flow, lb/ sec 

0.012832 

1.71468 


Table I. Minimum and Maximum Test Conditions 
of Test Matrix for Functional Relationship 

E = f(u,p,T„T w )\ Method I. 





Total Pressure, psi 

3.675 

29.4 

Static Pressure, psi 

3.44859 

26.2216 

Mach Number 

0.10 

0.95 

q , psi 

0.0066 

9.8543 

Total Temperature, °R 

540 

540 


457.433 

539.73 

Reynolds number/ foot 


itfffsfca 

Knudsen number 

0.00494 

0.10369 

Heater Power, Watts 

33.2 

2985 

Density 

0.01875 

0.13333 


0.013103 

1.1793 

| Overheat Ratio 

0.3 

1.0 


Table n. Mtnimnm and Maximum Test 
Conditions of Test Matrix for Functional 
Relationship E= f(M,Kn,a' w ); Method II. 
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Results an d Disc ussion 

Data obtained in the " Probe Calibration 
TunneT was used to calculate the sensitivities of 
several hot-wires using the two methods developed 
for CTA. The data were obtained such that the partial 
derivatives of the functional relationship such els 
E = and E = f{M,Kn,a‘ w ) were satisfied as 

dictated by the csdculus, i.e., while obtaining data by 
VEuying one independent variable, the other 
independent variables were held constEint ConstEmt 
Temperature Anemometers (CTA) were used in 
obtEiining the hot-wire data. Thus, the restriction of 
the calculus with respect to T w was satisfied 
automatically Euid the first functional relationship 
reduces to E = f(u,p,T t ) . Polynommls of appropriate 
degree were chosen to give the " best curve fit to the 
data to evaluate the partial derivatives required in the 
sensitivities equations. 

The test conditions in this study for 
E = f(u,p,T 0 ) case sire presented in Table I. Examples 
of curve fits to the data Eire shown in figures 6 (a-b) 
where the dependent vEiriable, measured voltage, is 
presented against each of the independent variable, 
while the other independent vEiriables were held 
constsint. Computing sensitivities in this case was 
simple and direct. The logarithmic derivative of the 
measured voltage to each independent variable gives 
the sensitivity of the hot-wire to corresponding 
variable as shown in equation (16). 

In the second set of data, a' w was chosen over 
T 0 to demonstrate that even though many wind 
tunnels do not have the capability of varying total 
temperature, the limitations of the tunnel can be 
overcome by using a' w . In obtaining the second set of 
data with respect to overheat ratio, all the hot-wires 
were operated at assigned overheat values. The test 
conditions for this case are given in Table II. The 
curve fit of measured voltage vs. vEirious independent 
VEiriables are presented in figure 7 (a-b). Here the 
sensitivities were computed using equations (24)-(26). 

To obtsun fluctuation quantities each hot-wire 
on the probe was assigned different overheat value in 
an attempt to make 5. * S p , a necessary condition for 
obtEiining fluctuations. However, no fluctuation 
quantities are presented in the present study. 


Velocity Sensitivity: 

As seen in figures 6-a and 7-a, the voltage 
increases rapidly at low velocities or Mach numbers 
and at transonic speed the slope of the curves tend to 
decrease. For the case E = f(M,Kn,a' w ) there was a 
tendency for dE/dM —* 0 at intermediate Mach 
numbers before increase as M -* 1. Thus, in 
correlating E vs. M at transonic speeds might result 
in zero or negative slope. This variation was eiIso 
observed in the data presented by other 

researchers®**®* * *. In order to correlate these data 
higher order polynomials were required. This could 
lead to large difference in the sensitivities at higher 
Mach number due to poor curve fitting at the end 
points. Precautions need to be taken in using the 
proper degree of polynomiEds because they could lead 
to unreasonable values for the sensitivities. One 
such difUculty in using higher order polynomial was 
presented in reference 9 where it lead to a hairpin 
like VEuiation of 5. with S p . 

Examples of the velocity sensitivity, S„, for 
both sets of data are presented in figures 8 (a-b) and 
generally, 5„ ranged from 0.05 to 0.40. In both cases 
S„ increased with increased u or M at low subsonic 
speeds and decreased at intermediate subsonic 
speeds. As S m approached sonic velocity, the 
sensitivities increased rapidly. 

A comparison of 5, for both the methods is 
presented in figure 9. It is expected that both 
methods, at identical test conditions, should give 
similar results. Most of the values for S m agreed 
within about 25 percent, however, some of the 
differences shown may be due to different degree of 
polynomial curve fits used in each technique. Also, it 
may be due to the differences in the operation of the 
Euiemometer for each method. In the first method, T. 
is held constant smd T c is varied. In the second 
method, T c is held constant and T w is vsiried. Due to 
this difference in the operation of the hot-wire the 
end-loss to the support could be different This 
difference in end loss could contribute to differences 
in all the three sensitivities. 


Density Sensitivity 

The measured voltages across the hot wires 
monotonically increased or decreased with increasing 
density or Knudsen number, respectively. At higher 
Knudsen number the slope chsuiged sign which was 
due to large change in Knudsen number and 
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insufficient data . A second degree polynomial was 
used to correlate the data. The density sensitivity, S p , 

varies from 0.1 to 0.4 and an example of S p for both 

the cases are presented in figures 10 (a-b). Except for 
the data where Knudsen number had a large change, 
the density sensitivities in both the cases were 
comparable with each other and a case is presented 
in figure 11. 

A comparison of velocity sensitivity, S„, with 
density sensitivity, S p , is presented in figures 12 (a-b) 
and shows that S u * S p for most of the data. However, 
few data points do approach S m = S p and S„> S p in 
some cases for E - f(M,Kn,a' w ) . The cases where 
S p > S „ indicates that the wire was in the slip flow 

regime. This observation is consistent with 
Spangenberg's data® as noted in reference 9. It was 
also noted that at higher Mach numbers and lower 
overheat ratio, a'„, S m approaches S p . As a' w was 
increased S, became greater than S p . 

Total Temperature Sensitivity 

For method II, the manufacturer's suggested 
values of a, and Morkovin 1 suggested values of >8, for 
tungsten wire were used to compute . a , was also 
computed using the data obtained during the test 
and this value was comparable with the 
manufacturer's provided value of a,. The data 
obtained in the present study was not adequate to 
evaluate fi,. 

An example of the data used to obtain the 
total temperature sensitivity, S ,., is presented in 
figures 6(b) and 7(b) and S z ranges from -0.2 to -1.0. 

However, the total temperature sensitivity obtained 
using the first method was always larger them the 
values computed using the second method. A few 
examples of these variations are presented in figures 
13 (a-b). This range of temperature sensitivity is 
consistent with the earlier studies®’^ from the 
authors. The temperature sensitivities using the 
overheat parameter may be more accurate because 
better curve fits were possible since more data points 
were available when using the E = f{«,p,T 0 ) method. 
This shows the possible advantage of using this 
method versus the other. 

A comparison of total temperature sensitivity, 
S t , between the two methods is presented in figure 
14. There is a significant difference between the 


values of for the two methods. This difference 
might be explained if the end losses are taken into 
account since they might be different due to the 
difference in the operation of the anemometer for 
each method. 

In general, the range of data available from 
the present study is greater than those obtained from 
previous studies®*®, particularly in regards to 
temperature information. In some of the earlier 
studies logarithmic functions were used to correlate 
the data which gave no variation of the sensitivity to 
its corresponding independent variable. Therefore, 
polynomial functions seem to be a better technique 
for correlating the data compared to logarithmic 
functions. However, careful analysis of the present 
data is required for a better understanding of the 
differences between the sensitivities, particularly the 
total temperature sensitivity, so that they could be 
applied for computing the fluctuation quantities. 

Conclusion 

Conventional methods for calibrating hot-wire 
probes use equations which include quantities that 
are difficult to evaluate. These quantities are wire 
length and /9, in the temperature - resistance 
equation for the wire material. This could lead to 
possible errors in obtaining the sensitivities of the 
heated wire to changes in velocity, density and total 
temperature. Methods were developed that either 
eliminated the need for the measurement of t and >8, , 
which are difficult to make accurately or reduced the 
effect of /?,. From a study of these methods the 
following conclusions can be made: 

1. The equation obtained for a CTA based on the 
assumption that E = f{u,p,T 0 ) resulted in an 
equation in which neither wire length or a, and 
/?, are required. This should result in a more 
accurate method for calibrating hot-wire probes 
in facilities where T c can be varied. 

2. The equation obtained for a CCA based on the 
assumption that E = f(u,p,T e ,T w ) results in an 
equation that requires the variation of four 
independent variables and is believed to be too 
complicated for routine use. The equation also 
requires the evaluation of /?, who's elimination 
was desired. 

3. The equation obtained for a CCA and CTA 
under the assumption that E = f(M,Kn,a' w ) and 
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T a = constant resulted in an equation where a, 
and Pi occurred in only one term in the 
equation for the total temperature sensitivity. 
Because of this the calibration of wires should 
be more accurate than those obtained using 
conventional methods. 

4. The sensitivities computed based on the 
methods presented using the data obtained in 
* Probe Calibration Tanner gives comparable 
results for S, and S p . However, there were 
significant differences in the values for . The 
differences in the sensitivities could be due to 
different degree of the polynomial curve fits 
used in each technique or due to the 
differences in the operation of the anemometer 
in each method. 

5. In general, the velocity sensitivity ranges from 

0.05 to 0.40, density sensitivity ranges from 
0.10 to 0.40 and total temperature sensitivity 
ranges from -0.2 to -1.0. 

6. Polynomial functions seem to be better for 
correlating hot-wire data compared to 
logarithmic functions. 
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Figure 1. Schematic Diagram of Heated Wire. 
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Figure 3. Schematic Diagram of Data Acquisition 
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Figure 4. Operational Envelope of the " Probe 
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Abstract 

Much has been written about the improbability and impossibility of obtaining solutions to the mean 
square equation for constant current anemometiy in subsonic slip and transonic flows. For these flow 
conditions, the fluctuating voltage across a wire mounted normal to the flow is a function of velocity, density, 
and total temperature. In principal, the fluctuations of these quantities can be measured, however, to date there 
are no known acceptable solutions to the mean square equation in these flow regimes. In this study, data 
presented in the 1950's by Spangenberg were used to show that, for a large region in the Nusselt number - 
Reynolds number or Nusselt number - Knudsen number regimes, there exists the possibility of obtaining 
solutions to the mean square hot-wire equation. These data were used to compute the sensitivities of the heated 
wire to changes in velocity, density, and total temperature; indicate regions where the velocity and density 
sensitivities were different (a condition required for a solution to the mean square equation); and show a second, 
necessary condition for a solution. Examples of fluctuation and mode diagrams for subsonic compressible and 
subsonic slip flows are also presented under the assumption that there are solutions to the hot-wire 
anemometry equation. 

Nomenclature 

a 0 _a. coefficients in equations (22) - (27) 
b 0 -b m coefficients in equations (24) - (26) 
c 0 - c. coefficients in equation (26) 
a,b,c &.d constants in equation (A-4) 

A quantity given by equation (20) 
c f specific heat at constant pressure 

c, specific heat at constant volume 
d wire diameter 

e' fluctuating voltage across the sensor 
E w mean voltage across sensor 
f,g,h constants in equation (A-4) 
h coefficient of heat transfer 
/. current through the wire 
k c thermal conductivity of air based on T„ = 0°C 

K <?logr./<3ogf^ 

Kn Knudsen number 
l wire length 
L characteristic length 
m mean mass flow 


M 

Nu 

n 

P, 

Q, 
1 

K 


S 

s. 

s 


T 

r 


Mach number 
Nusselt number 
number of data points 
electrical power to wire 
heat transfer from wire 
sensitivity ratio = S./Sj- 
resistance of wire 

wire resistance at reference temperature 

sensitivity ratio = /S ? 

sensitivity ratio = S p fS Xr 

mass sensitivity 

velocity sensitivity 

density sensitivity 

total temperature sensitivity 

temperature 
total temperature 
velocity 

general dependent variable 


This paper is declared a work of the U.S. Government and is 
not subject to copyright protection in the United States.. 
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a 

a, 

P 

P\ 

A 

n 

e 

e 

r 

p 


general independent variables 
-(^og/./^log«.) 

linear temperature - resistance coefficient 

of wire 

a(r-\)M 2 

second degree temperature - resistance 
coefficient of wire 
mean free path 

recovery temperature ratio, T^/T, 
temperature parameter, T w /T 0 
angle between probe axis and plane sound 
wave 

specific heat ratio, c r jc r 
density 

temperature loading parameter, (T". - 7^J/7^ 
normalized fluctuation voltage ratio = ^e'/ 


Superscript 

indicates instantaneous value 
RMS value 
mean value 


Subscript 

adw adiabatic wall condition 

w wire condition 

oo free stream or static condition 

o based on T = 0 °C 

ref reference condition at T 0 = 0°C 

to vorticity 

er entropy 

sf far field sound 

Introduction 

Constant current anemometry (CCA) has been 
extensively used to measure fluctuations in 
incompressible and compressible flows In 

compressible flow the heated wire is, in general, 
sensitive to changes in velocity, density, and total 
temperature and fluctuations in these quantities can, 
in principal, be measured'*. In supersonic flow it has 
been found experimentally that the mean voltage 
measured across a heated wire is only a function of 
changes in mass flow and total temperature®. Using 
this result for supersonic flow, Kovasznay developed a 
graphical technique for determining fluctuations of 
mass flow, total temperature, and their correlation®. 


Very few attempts have been made to measure 
fluctuations in subsonic compressible flow using a 
CCA because of the complexity that exists when the 
heat transfer from the wire is a function of three 
independent variables. 

Some efforts have been made to measure 
fluctuations in subsonic compressible flow using a 
constant temperature anemometer (CTA). For these 
investigations a three-wire probe was used and the 
voltages across the wires were digitized and a system 
of three equations solved to obtain the instantaneous 
changes in velocity, density, and total temperature as 
a function of time. Statistical techniques were then 
used to obtain fluctuations and correlations of 
interest^. 

For compressible flow; where the wire is 
sensitive to velocity, density, and total temperature; 
equations for CCA are usually derived using mean 
square values. This results in a single equation with 
six unknowns. Much has been written about the 
improbability'* and impossibility** of obtaining 
solutions to this mean square equation. Solutions to 
this equation requires that the velocity and density 
sensitivities of the heated wire be sufficiently different 
In this report, data published in the 1950’s by 
Spangenberg 9 was used to show that there exist a 
large region in the Nusselt number- Reynolds number 
or Nusselt number-Knudsen number regimes where 
the velocity and density sensitivities appeared to be 
sufficiently different to permit suitable solutions to the 
mean square equation. It was further shown that 
solutions were possible only when the velocity and 
density sensitivities were non-linearly related in 
specific manners. Again Spangenberg* s data were 
used to show that this second condition can often be 
satisfied, further indicating that solutions to the mean 
square equation are possible. Examples of 
fluctuation and mode diagrams for subsonic 
compressible and subsonic slip flows are also 
presented under the assumption that there are 
solutions to the hot-wire anemometry equation. 


Spangenberg** Data 

Spangenberg’s data were presented in terms 
of = f{M,p,T w ) for T c = constant. Data were 

presented for several wire diameters and wire lengths. 
The data to be considered herein were obtained using 
a FT- 10% Rh wire having a diameter of 0.00015 inch 
and a length of 2.80 mm. Data were limited to the 
following conditions: 
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0.05 < M < 0.95 
0.0004 < p, gm/cm 3 < 0.0012 
0.22 1< r. <2.192 

The Nusselt numbers presented were 
corrected for heat loss to the supports of the wire. 
This corrected Nusselt number, , is suitable for 
comparing heat transfer results from different 
experiments using different wire materials and 
different t/d's. For computing fluctuations, the 
sensitivities of the wire must be obtained using the 
uncorrected Nusselt number, Nu For the results 
presented herein, A values were converted to Nu e 
values using the method presented by Spangenberg. 

Examples of Spangenberg' s data are 
presented in figure 1. In figure la the uncorrected 
Nusselt number is presented as a function of Mach 
number for various values of density for a given value 
of r„. The Nusselt number increases rapidly with 
Mach number at veiy low values. The slopes of the 
curves decrease somewhat as the Mach number 
increases. At the higher Mach numbers (0.6 < M < 
0.95) the Nusselt number often reaches a peak prior 
to decreasing. This variation of Nu o with M at the 
higher Mach numbers often results in zero or negative 
slopes for Nu 0 vs M. The Mach numbers at which the 
slopes of the curves reach zero increase with 
increasing density. This variation of Nu 0 vs M was 
also observed in data presented in references 10 and 
11. Polynomials were used to fit curves to the data 
from which the required derivatives were obtained, 
and a fourth degree equation was required to obtain 
"good" fits. The variation of Nu t with p , shown in 
figure lb, is monotonically increasing with increasing 
density. A second degree equation was used to fit 
curves to the Nu 0 vs p data. 

The variation of Nusselt number with respect 
to r„ is presented in figure lc and this variation of 
Nu Q is rather complex and nonlinear. The Nusselt 
number can often increase then decrease with 
increasing values of r.. A closer look at Nu a vs r» 
data revealed that, for most of the data, the values of 
Nu 0 decreased from 0.629 Sr,S 1.629 but most often 
increased or tended to increase from 1.629 s r„ S2.192 . 
This variation of Nu t over this later interval of r„ was 
present in about 80 percent of the data and had a 
profound affect on the relationship between the 
velocity and density sensitivities. It is not known if 
this variation can be supported by the limited amount 
of data which is available. If third degree curves were 
fitted to the data, there was a significant change in the 


derivatives for most of the data at the higher values of 
r„. Because of these changes, the values of Nu t were 
fitted to values of r. using both second and third 
degree equations. The results obtained using these 
two curve fits on the variation of the velocity and 
density sensitivities will be presented subsequently. 

Data presented by Baldwin in Ref. 10 also 
showed a similar variation of the Nusselt number with 
r„. This complex and nonlinear variation of Nusselt 
number with r„ or overheat is probably the reason 
that solutions to the mean square equation for CCA 
might exist 

In general, the scatter in the values of N u 0 
with respect to M and p are acceptably small; the 
apparent scatter of Nu 0 with respect to r, is 
somewhat larger. This is probably due to the 
relatively small change in Nu c with respect to r„ or the 
limited number of values of r„ that are available. 


Theoretical Considerations 


Because of the manner in which Spangenberg 
presented his data, equations must be derived for the 
change in voltage across the wire in terms of the 
sensitivity of the wire to changes in velocity, density, 
and total temperature. The equation for the heat 
transfer from the heated wire is generally given as: 

2= ndthfo-Tj) (1) 

For steady state conditions, power supplied to the 
heated wire equals the aerodynamic heat transfer 
from the wire. Therefore: 


K = iIk = *tk c (T„ - r,T o )Nu o (2) 

Making the small perturbation assumption, the 
change in R, can be expressed as: 

(l - 2e)d log R, = c/log Nu o + d\o% r, + d log T„ (3) 


where: 


£= a°g/„ 

dog 

For Spangenberg's data: 
Nu e = f(M,p,T.) 


( 4 ) 


( 5 ) 
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Note that there is a direct relationship between 
density and Knudsen number. Therefore, the 

Knudsen number could be used as independent 
variable. The change in Nu c is: 

d\o%Nu 0 - ogr. («) 

mogp M (TiQ% T w 

It is desired to obtain the final equation in the 
following form: 

E. = f(u,pjJ (7) 

The change in M is: 

</log M = -{dlogu- i</log T„ j ( 8 ) 

and 

dlog r. = — JlogT. - —dlog rj~—d\og f 0 (9) 
r. r„ r. 

The recovery temperature ratio can be written as: 

»7= f(pM) (10) 

and the change in T) is: 

d\ogrj= ^ log p + — jf ^ log u - log 7 1 (11) 

^logp 2 y 


Morkovin^ gives the following equation for the 
changes in E w due to changes in u , p and T, for a 
CCA: 


«Mog£„ = -S a d\ogu- S /> d\o%p+Sj'd\ogT <> (16) 


For the present case the sensitivities can be obtained 
using equations (3) - (15) and are given by: 


5. 



V dlog tj ( dog Nu a i j'j 

dlog M 

r„ SogM[ dlogr w J 

dlog Nu. 

t) dlogrjf dlog Nu 0 ( j'j 

dlogp 

r„ dlogp\ dlog r„ J 


(17) 

(18) 


S, = j- J. SsslV u k + + 1 Y JL. ^ v. _ i)_ ] (19) 

1 [ 2a <?logr. J^2a<?logAS J ^og r» J 

where 



For Spangenberg’s data tj varied with M and 
plots of >7 vs M were presented in his report The 
values for a, and /?, for the wire material were also 
presented. Therefore, after the partial derivatives in 
equation (17) - (20) were determined the sensitivities 
to changes in u , p , and T o could be calculated. 


The temperature of the wire must be obtained from 
the measured resistance using: 

~ = 1 + «.( r - ' T J + Afa - T„ f ? (12) 

K *<s 

Therefore: 

T. = AK) (13) 

and 


The value of e cannot be evaluated from 
Spangenberg's data because the voltage of the battery 
used is unknown and there are several variable 
resistors, having unknown resistances, in the 
electrical circuit. Because of these problems the 
absolute values of the sensitivities could not be 
obtained. The values of e does not, however, 
influence the ratio of the sensitivities, the important 
parameters, since the quantity A cancels out when 
the ratios are formed. 


dlogT, = KdlogR . 


(14) 


The Partial Derivatives 


Finally, Ohm's law can be used to express the wire Spangenberg's data were measured and 

resistance in terms of the voltage across the wire as tabulated in a form that is necessary to obtain the 
follows: partial derivatives of the dependent variable with 

respect to one independent variable while holding the 
» n _ d\ogE w remaining independent variables constant. These are 

(1 + e) the only data known to the authors where 

measurements were obtained in this necessary 
manner. For the present report, polynomials were 


4 



used to fit the data to obtain the desired derivatives of 
Nu c with respect to the independent variables. The 
degree of the polynomial was chosen to give the 'best 
fir to the data. In general, this method gives: 

V=f{x,y,z) (21) 


Tti» M^m Sq uare Equation 

In the past hot wire data, using a CCA, were 
obtained from mean square values. Following 
Kovasznay, squaring equation (16), dividing by S'?, 
and taking the mean gives: 


where the function was assumed to be: 
V(x,y,z) = a„ + a t x + a 2 x 2 + ■■■a,x’ 


( 22 ) 


# 2 = 



(29) 


and: 

a m = f(y,z) m = 0,1,2,. ..n (23) 

Using polynomials, the values for the a 's can be 
expressed as: 

a m - K + V + + - b y (24) 

where 

b m =f(z)-, «= 0,1,2 n (25) 

In a similar manner: 

*» = c o + c i* + Cj2 J +••■ V" (26) 

From equation (22), the derivatives becomes: 

dlogV = a t x + 2 a 2 x 2 + ...na.x" 

<^og x a 0 + a, x + a 2 x 2 + ...ape" 

The Sensitivities 



This is a single equation in six unknowns, 


u p j u'fi u'T: fiT; 
u’p’T , ;* up'uT/ pT,; 


2 

Since ^ , q and s are 


functions of the overheat of the wire, one can, in 
principal, operate a single wire at six overheats and 
solve a system of six equations for the fluctuations 
and their correlations. The unsuccessful attempts 
that were made in the past to obtain solutions in this 
manner were thought to be due to the mean flow 
calibrations not being accurate enough or the 
sensitivities not being sufficiently different to permit 
accurate solutions 4 . Demetriades® stated that no 
solutions to equation (29) are possible unless the 
independent variable occurs to at least the fifth 
power. It is shown in the Appendix that if s is 
assumed to be a function of q , solutions to equation 
(29) are possible if, after substituting the functional 
relationship for s into the equation, the resulting 
equation has six terms. This result requires that 
values for s have very restrictive, non-linear 


variations with respect to q as the overheat of the 
wire is changed. Table I presents conditions where 
solutions to equation (29) exist based on s being a 


simple, finite, power series of q . 


After the partial deviates were determined, the 
sensitivities were obtained using equations (17) - (20). 
In general: 

S.-S^ = /(M,aO (28) 

For these calculations it was assumed that e=0. 
Plots of the sensitivities will not be presented since it 
will be shown subsequently that the ratios S u / S z and 
S"/Sj. are the important parameters. At some values 
of the Mach number, S m = 0 and the wire, at these 
conditions, is not sensitive to changes in velocity. 


Figure 2 is presented to show to what extent 
q * s for r. = 2.192. From the figure it can be seen that 
s is, in general, greater than q and the difference 
can be large since q can often approach zero at the 
higher values of M. At the higher values of M the 
sign of q can change and negative values are not 
presented in the figure. 

For Spangenbeig's data, which is considered 
herein, the Mach number ranged from 0.05 to 0.95. 
The equation given by Morkovin for the Knudsen 
number is: 


7.75x10 9 
fd 


(30) 
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where p is in gm/ctn 3 and d is in cm . Using this 
equation, the value for the Knudsen number for these 
data ranged from 0.017 to 0.051. Reference 12 
suggest that continuum flow occurs for Kn < 0.001 and 
slip flow occurs for 0.001 <. Kn £ 2 . Other references 
suggest that continuum flow occur for Kh <0.01. 
These quoted boundaries are not sharply defined, but 
if it is assumed that slip flow occurs for Kn > 0.01, all 
of Spangenberg’s data considered herein is in the slip 
flow regime. This result is in agreement with the data 
presented in figure 2 where SJS^ * SJSj. , however, a 
few data points do approach S u /Sj. = S p fS^ . 


Solution to Mean Square Equation 

Solutions to equation (29) can be obtained by 
using the method of least squares. Applying this 
method results in the following matrix equation: 


£* 4 Zq 2 ? 

Zq 2 

w 

Sf 3 

Zq 2 i 

Zs* 

zP- 

ZqP 

I-P 2 

zP 


n 

Zqs 

29 

Zs 



Zq 2 P 

29 2 r 

ZqP 




Zq 2 

Z‘q 





zP 


id Ip ) 1 
(W? 
2 0JL 

« p 

T 0 u 

£IL 

PT. 


\-2 


zsE 

2fV| 

= 

ZqjP 

Ztf 2 


( 31 ) 
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In equation (31) q , a and ^ are functions of the 
overheat parameter, r„. Solutions to equation (31) 
can be obtained as follows. Obtain a set of values for 
~2 

q , a and 4> for a suitable number of values for r„, 

~2 

say 10 or 12. Fit suitable curves to q and 4> vs r„. 
Also obtain curve fits for a = f{q) that satisfies the 

requirements noted in the Appendix to determine if a 
solution to equation (31) is possible. The * best curve 
fit to these data using the acceptable coefficients can 
be used to relate a to q and ultimately a and r„. 

T 

Using the equations obtained for q , a and 4> vs r w , 
compute a set of values, say 20, to be used in 
obtain the fluctuations and their 

correlations. 


An additional datum point can be obtained for 
s and q by using the limiting values for these 
variables. These equations are: 


<3og>7 


1 <5\ ogM 

of 1 dog? jV ' 

(32) 

\la c\ogM ) 


ologr? 


— aaa— r;r..o 
i <aog>? ^ 

(33) 


1 flog 17 "j’j 

cftogM J 


An indication of the accuracy and suitability of these 
results can be obtained by computing the condition 
number f or the 6x6 matrix in equation (31). The 
smaller the condition number the more accuracy can 
be expected from the calculated results. 

No fluctuating voltages were presented by 
Spangenbetg, therefore, no fluctuating quantities 
could be calculated. His data were used to obtain 
relationships between q and a , to calculate the 
elements of the matrix in equation (31), and to obtain 
the condition numbers for these matrices. This 
procedure serves to illustrate that solutions to 
equation (31) are possible and give some indication of 
the accuracy of possible solutions. 

Examples of the variation of s with q 
obtained from Spangenberg's data is presented in 
figure 3. The case where the Nusselt number was 
assumed to be a cubic relationship with respect to r. 
most often resulted in a hairpin like variation of a 
with q . This type of variation does not appear 
reasonable, therefore, a curve is also presented in the 
figure where Nu„ was assumed to be given by a 
second degree equation. This assumption resulted in 
what appears to be a more reasonable variation of / 
with q . It must be again noted that the behavior 
obtained with the cubic curve fit might be due to the 
limited amount of data available. 

There were several sets of data where Mi, 
could be fitted with a cubic equation with respect to 
r„, without the derivative of Nu e increasing at the 
higher values of r„. Examples of a vs q for these 
cases are presented in figure 4. The powers of q 
presented in Table I were fitted to the a- q data, and 
the powers that gave the ’ best curve fits!' to the data 
are presented in the Table n. Most of the curves fitted 
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to the data appear to be good based on the F? values 
except for the case M = 0.60 . Curves fitted to the data 
using the values of the powers presented in Table II 
are given in figure 4. The equations with these 
powers were used to relate a to q and in turn to r. 
Twenty values for q were selected, the values for a 
computed, and the elements for the 6x6 matrix in 
equation (31) calculated. Using the elements of the 
matrix, the condition numbers for the cases presented 
in figure 4 were determined and these results are 
presented in the Table II. These large numbers 
indicate that the 6x6 matrix in equation (31) is veiy ill 
conditioned. It is not known if these results are 
general or limited to the cases investigated which has 
a limited number of overheats. There are techniques 
that can be used to improve solutions for equations 
having ill conditioned matrices It is not known if 
these techniques cam be used to assure suitable 
solutions to the hot wire equation. There are cases 
where suitable solutions can be obtained to sets of 
equations where the matrix results in large condition 
numbers^. It must be demonstrated that 
satisfactory results can be obtained for the present 
equation. Thus, the necessary condition for a solution 
to the mean square equation has been established; it 
remains to be demonstrated that sufficient conditions 
can be established for suitable solutions. 


Fluctuation and Mode D<«yram« 

Fluctuation Diagram! 

It is recognized that, in general, if solutions to 
equation (29) can be obtained using the method of 
least squares then the fluctuations and their 
correlations can be obtained without considering the 
fluctuation diagram. However, the fluctuation 
diagram can often be used to advantage in 
determining the dominant mode in a flow^. 


Before considering the fluctuation diagram for 
subsonic compressible flows, consider the case for 
supersonic flow where the velocity and density 
sensitivities are equal. Under these conditions 
equation (16) becomes: 



(34) 


Dividing the above equation by S z , squaring, and 
taking the mean gives: 



There is an additional problem that could 
have an adverse effect on solutions to equation (29). 
Equations (32) and (33) show that a and q , in 
general, cannot be zero at the same time. This means 
that conditions required to obtain the total 
temperature fluctuations do not lie on the curve 
defined hy s= /(<?). Therefore, the total temperature 
fluctuations must be obtained by extrapolation to the 
condition where q= s= 0. It must be noted that this 
condition is not determined by the accuracy of the 
mean flow calibration data but by the relationships 
used to obtain equation (29). This problem is not 
limited to the subsonic compressible flow case, but 
could also exist in supersonic flow if dlogrj/d\ogRe * 0. 
In supersonic flow the extrapolation is in two 
dimensional space, whereas in subsonic flow the 
extrapolation must be in three dimensional space (see 
Fluctuation and Mode diagram section). 
Extrapolation, in three dimensional space, of a 
nonlinear function could result in highly inaccurate 
solutions. At the present time it is not known whether 
other relationship could be used to obtain values for 
a and q such that this problem can be rectified. 
Possible remedy might be found in the techniques 
presented in reference 11. 


This is a single equation with three unknowns 
ih f fti * V 1 

».e., — , I. In principle, a single wire can be 

m T 0 mT 0 J 

operated at three overheats and a system of three 
equations solved for the fluctuations and their 
correlation. Kovasznay^ noted that this technique 
would probably result in inaccurate answers. He 
suggested that the wire be operated at several 
overheats and a graphical technique be used to obtain 
the fluctuations. Equation (35) is an equation of a 
hyperbola where the asymptote gives the mass flow 
fluctuation and the intercept of the curve with the i - 
axis gives the total temperature fluctuation. Note that 
the method of least squares can also be used to obtain 
solutions to equation (35). An example of a general 
fluctuation diagram for supersonic flow is presented 
in figure 5. 

Kovasznay^ demonstrated that fluctuations in 
compressible flow are composed of three basic 
fluctuations which are vorticity, entropy, and sound. 
He denoted these three basic fluctuations as modes. 
Vorticity is a velocity fluctuation that is usually 
described as turbulence and has no pressure or static 
temperature fluctuation. Entropy consists of static 
temperature fluctuation or temperature spottiness. 
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The sound field consists of isentropic fluctuations 
associated with pressure fluctuations and include 
density, velocity, and temperature fluctuations. If 
fluctuations are considered to be composed of a single 
mode then the fluctuation diagram was noted as a 
mode diagram. 


Next, assume that the only fluctuations are 
entropy. This results in: 



If it is assumed that the velocity and density 
sensitivities are sufficiently different and sufficiently 
non-linear so that solutions to equation (29) are 
possible, what are the characteristics of the 
fluctuation diagrams in subsonic compressible flow? 

.2 

In equation (29), l is a function of q and s , 
therefore, the fluctuation diagram exists on a three- 
dimensional surface, a hyperboloid, rather than a 
plane as in the case when 5. = S p . The locus of points 

of the fluctuation diagram on the surface of the 
hyperboloid will depend on the relative changes in q 
and a as the overheat of the wire is changed. 
However, the important information exists in the ~<f>-q 
and planes. For example when j= 0, equation 
(29) reduces to an equation for a hyperbola in the 
j>-q plane where the asymptote gives the velocity 
fluctuations. If q= 0, again equation (29) reduces to 
an equation for a hyperbola in the j-a plane and the 
asymptote represents the density fluctuations. When 
q and a are zero, the intercept on the ~<f> -axis gives 
the fluctuations for the total temperature. In planes 
parallel to the q-a plane, the locus of points of the 
fluctuation diagram is governed by the velocity and 
density fluctuations and their correlation. The cross 
product term, qa, requires a rotation of the axis before 
the characteristics of this locus can be identified. 

Although the fluctuation diagram exists on 
the surface of a hyperboloid, the fluctuations can be 
determined from the intersection of the hyperboloid 
with the } \-q and planes. Because of this, the 
fluctuation and mode diagrams will be defined as the 
traces of these intersections in the noted planes. A 
general schematic representation of the fluctuation 
diagram for equation (29) is presented in figure 6. 


Finally, assume that the only fluctuations are 
far-field sound. Then: 


«*_ 1 ( P’ 


fill toss ; (38, 

l p )<r P y\P) 4 r - 7 \P)^ 

Ht)/ f),‘ 39 > 


A schematic of the mode diagrams for 
supersonic flow is presented in figure 7 for vorticity, 
entropy, and far-field sound. In all cases the mode 
diagrams are degenerate hyperbolas. This figure is 
presented in order to compare the mode diagram in 
supersonic flows with those in subsonic flows. 

If it is assumed that in subsonic flow the 
fluctuations are only due to vorticity, equation (29) 
gives: 


^• 2 = {q-pf{~^ (40) 

Results obtained using this equation for the case 
where (w'/u)^ = 0.01 are presented in figure 8. The 
mode diagram for vorticity is identical to the one for 
supersonic flows except that S m is replaced by 5. and 
the important part of the mode diagram lies in the 
~<p-q plane since p/p=0. Equation (40) shows that 
when $ = 0 , then q = p. Also, when q = 0 then: 
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which represents the total temperature fluctuation 
due to vorticity. 


Mode Diagrams 

First, assume that the only fluctuations are 
vorticity. Then: 



Next assume that the fluctuations are all 
entropy. Then equation (29) becomes: 
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This equation is identical to the one for supersonic 
flow except that S m is replaced by S p and the mode 

diagram lies in the plane. 

An example of the mode diagram for entropy 
when {f m /T m ) =0.01 is presented in figure 9. Here 

again, only the +-s plane is presented since uju = 0 . 
When ^= 0 then 4 = -a, and when 4 = 0 then: 


<(>= a 


f f \ 


(43) 


which represent the total temperature fluctuation due 
to entropy fluctuations. 

Now assume that all the fluctuations are far- 
fleld sound (i.e., X/L« 1). The equation for this case 
is: 


2 

<*' 2 = y^ki 1 (~p-) ^"^ COs6 ' + M r_ ( 44 ) 

Equation (44) represents a three-dimensional surface 
in the f , q , s coordinate system and the mode 
diagram, in general, lies on this surface. The 
intersections of the surface with the }>-q and ~4>-s 
planes are defined, herein, as the mode diagrams. 
When 4=0 the above equation represents the mode 
diagram in the q plane and when 9 = 0 the mode 
diagram is in the ^-4 plane. An example of the 
sound mode is presented in figure 10 , where the }>-s 
plane has been rotated into the plane of the paper. 
The sound mode is the only mode that has non-zero 
values in the two planes, j-q and since sound 
has both velocity and density fluctuations. All of the 
mode diagrams described above are degenerate 
hyperbolas in the ~j-q and j-s planes. 


• It was shown that the mean square equation for 
CCA has possible solutions provided that the 
functional relationship which adequately relates 
4 to q when substituted into the equation results 
in only six terms. Again Spangenberg’s data were 
used to show that this condition can often be 
satisfied. 

• The application of the method of least squares to 
the mean square equation resulted in a 6 x 6 
matrix that was ill-conditioned for the data used. 
Methods to circumvent this problem, if it is 
general, must be found before reliable solution 
can be obtained. 

• If a solution exists, the fluctuation diagram was 
shown to exists on a three dimensional surface - a 
hyperboloid. 

• The intersection of the hyperboloid with two of the 
co-ordinate planes are hyperbolas, where in one 
plane, (9 - ^) , the asymptote gives the velocity 

fluctuation and in the other plane, (4 - ^) , the 

asymptote gives the density fluctuation. The 
intersection of the hyperboloid with the i axis 
gives the total temperature fluctuation. 

• The mode diagrams for vorticity, entropy, and 
sound are degenerate hyperbolas in two of the co- 
ordinate planes. 
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Concluding Remarks 

From the present study of hot-wire 
anemometxy for a constant-current anemometer 
applicable to subsonic compressible and subsonic slip 
flows, the following concluding remarks can be made: 

• Based on Spangenberg’s data, there appears to be 
a wide range of M and Kn available where a 
solution to the mean square equation for constant 
current hot-wire anemometiy in subsonic flow 
might be possible based on the fact that S y *S p . 
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Appendix 


A Note on the Mean Square Equation for CCA 


The necessary condition required to assure a 
solution to the hot-wire equation for the constant 
current anemometer can be developed as follows. The 
equation in terms of mean square values is: 




where 



x S r 
and s=~ 


(A-2) 


In equation (A-l) the quantities l , q , and s 
are functions of r„. Therefore s can be considered to 
be function of q . Assume that solutions are sought 
using the method of least squares and j = f(q); 

substitute this functional relationship into equation 
(A-l): 


11. Nagabushana, K.A., and Stainback, P.C., 1994, "A 
Rational Technique for Calibrating Hot-Wire Probes 
in Subsonic to Supersonic Speeds," AIAA-94-2536, 
A1AA 19th Aerospace Ground Test Conference, 
Colorado Springs, Colorado. 

12. Stadler, J.R., Goodwin, G., and Creager, M.O., 
1951, "A Comparison of Theory and Experiment for 
High-Speed Free-Molecule Flow," NACA Rep. 1032 
(Supersedes NACA-TN-2244). 

13. Strang, G., 1976, "Linear Algebra and its 

Application,’ Academic Press. 

14. Ojalvo, I.U., 1991, "Improved Solution for Ill- 
Conditioned Algebraic Equations by Epsilon 
Decomposition," AIAA Journal, Vol. 29, No. 12, pp. 
2274-2277. 


(a ~ 3) 

’ Any functional relationship which adequately 
relates q and s and results in equation (A-3) having 
six terms will satisfy equation (A-l) and result in a 
solution. This can be illustrated by assuming the 
following equation for s : 

s= a + bq * + cq* + dq k (A-4) 

Substituting equation (A-4) into (A-l) and collecting 
terms gives: 
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Table I 



( AS ) 


Possible Solutions to Mean Square Equation for 
CCA: 


* 2 ' - Mif (f f 


•- /(?)“ a* bq* + eg* + dtf 


A solution to this equation is possible for any 
values of /, g, and h provided that the resultant 
equation, after the substitution of values for /, g, 
and h into the equation, has 6 terms. Some 
examples of values for / , g , and h which results in 
an equation which has six terms are presented in 
Table I. However, when equation (32) is considered, 
the non-integer values for / cannot be used to 
determine s in the vicinity of r w = 0. This is due to 
the limiting value of q being negative. This results in 
the values of / being restricted to the values 
presented in the first two rows of Table I. The major 
problem in obtaining a solution to equation (A-l) is 
evaluating the powers of q in equation (A-4) which 
will give the best fit to the mean flow calibration data 
relating s to q . This can be done by using standard 
statistical techniques. 


/ 

g 

* 

Terms in q 


/ 

i 

0 


f = 3,4,5,... 

/ = -2,-3, -4,... 

/ 

0 

0 


/ = 3,4,5,... 

/ = -2, -3, -4,... 


1 

0 


/ a non - ml *ger value! 
!</ 

' 

0 

0 


/ ® non-mtffer vokte* 
1 </ 

/ 

1 

0 

fWVw’w 1 * 0 

0</<l 
f* 0.5 

L 

0 

0 

WWm'hW 

0</<l 
f* 0.5 


TABLE II 


Possible Solutions To Mean Square Equation For 
CCA Based On Boangenbexg's Data: 


~~W~ 

P 

/ 

g 

h 

R 2 

Condition 

Number 

0.05 

0.0012 

3.0 

T" 

~0~ 

0.9936 

1.97 x 10 7 

O.iO 

0.0012 

3.3 

i 

0 

0.9966 

2.92 x 10 8 

0.20 

0.0012 

3.0 

i 

0 

0.9933 

1.92 x 10 8 

0.40 

0.0012 

3.0 

i 

0 

0.9997 

1.09 x 10 7 

1 0.60 

0.0012 

3.0 

i 

0 

0.9718 

6.54 x 10 8 
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ADDENDUM 

The present analysis of constant current 
anemometry in subsonic slip and transonic flows 
using the least square technique appeared to show 
that, based on Spangenberg's data, measured 
fluctuations would probably have substantial 
errors based on the large condition numbers 
calculated for the 6x6 sensitivity matrix. This was 
true even though the necessaiy conditions required 
for a solution to the mean square equation 
appeared to be satisfied. Because of this, 
additional studies of the problem were undertaken 
in attempts to determine the possible reason for 
the large condition numbers and to search for 
possible methods for reducing their magnitudes. 

From an analytical point of view, consider 
equation (29). One serious problem which could 
affect the ability to obtain accurate solutions to this 
equation can be noted by considering the variation 
of the sensitivity ratios, q and s , with respect to 
changes in the overheat of the heated wire and 
their limiting values which are give in equations 
(32) and (33). 

In subsonic slip and transonic flows, the 
values of s are, in general, larger than those for q . 
Both of these quantities approach small limiting 
values as the overheat of the wire approaches zero. 
However, except for very low densities, the limiting 
values for s are smaller in magnitude than those 
for q as the overheat approaches zero. Because of 
this, there is a range of overheats in the vicinity of 
zero where q and s are approximately equal and 
the locus of points for q and s in the q-s plane 
crosses the q = s line. Because of this, the mean 
square equation reduces from the case of a single 
equation with six unknowns to a single equation 
with three unknowns. In other words as the 


overheat is reduced, the mean square equation 
becomes degenerate from a three dimensional 
problem to a two dimensional one and this could 
result in an inability of the equation to give 
accurate results. Therefore, the present method 
would probably result in answers having large 
errors and accurate results cannot be expected 
unless some other suitable technique can be found 
to solve the mean square equation. Note that this 
problem in the limit is due to the heat transfer 
equation used and not due to any errors in the 
mean flow data. Is the heat transfer equation used 
adequate? If it is not, can another technique be 
found to circumvent the problem of obtaining 
accurate answers to the mean square equation? 

After the wire is calibrated and the 
sensitivities obtained and the RMS and mean 
voltage measured across the wire, the search for a 
possible solution to the mean square equation 
becomes one of geometry. The trace of the values 
for q and s in the q-s plane and the 
corresponding values of 4 > must be used to define 
the hyperboloid from which the fluctuations and 
their correlations can be accurately determined 
since the fluctuation diagram lies in its surface. In 
general, q and s exist in only one quadrant of the 
q-s plane as the overheat of the wire is changed, 
and the hyperboloid must be determined from this 
limited amount of data. Table II shows that 
attempts to do this resulted in matrices with large 
condition numbers indicating that any answers 
obtained would be subject to large errors. These 
errors might be reduced if additional information 
could be obtained. 

Is there a reasonable method that can be 
developed to supply additional information to help 
define more accurately the proper hyperboloid 
thereby producing more accurate solutions to the 
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mean square equation? If the functional 
relationship between q and s are valid in one 
quadrant, values for q and s can be obtained in a 
second quadrant by extending values from the 
equation into the second quadrant This technique 
reduces the condition numbers of the matrix by 
about two orders of magnitude. There are, 
however, no corresponding values for ^ in this 
second quadrant A method to circumvent this 
problem is to assume ^ is a function of q in the 
first quadrant then extend this relationship into the 
second quadrant Using this technique the 
condition numbers of the sensitivity matrix were 
reduced to the values presented in Table II. These 
condition numbers indicate that significant 
improvements in possible solutions can be 
obtained, however, the condition numbers are still 
rather laige. Because of the symmetry of the 
hyperboloid, the above technique might be 
acceptable. However, it should be noted that the 
coordinates of the hyperboloid is not symmetric 
with respect to the q , 3 , <f> coordinate axes. 

Based on the present paper and the 
Addendum, it appears that the possibility exist for 
obtaining solutions to the mean square equation 
for CCA in the subsonic slip and transonic flows. 
The accuracy of the results are , however, subject 
to question. In any case, the possibility of a 
solution suggest that data, including fluctuation 
measurements, should be obtained in an attempt 
to further investigate the problem. 


ERRATA 

In general, the values for q and s are 
negative and equation (16) is written so that the 
terms in equation (29) are all positive. For the 
present report the negative values for q and s 
were transformed to positive values. However, 
most of the equations were presented in the non- 
transformed coordinate system. This fact was not 
noted in the report and could lead to some 
confusion. Equations (31) - (33) are written in the 
non- transformed system. 

Page 1 1 of the Appendix 

However, when equation (32) is considered, 
the non-integer values for / cannot be used to 
determine s in the vicinity of r. equal zero. This is 
due to the limiting value of q being negative. 

The above should read: However, non- 
integer values cannot be used since the fluctuation 
diagram exist on the surface of a hyperboloid and 
this surface must be defined in all quadrants. 

The value 3.3 in Table II should be 3.0. 

The notations for the coordinate axes in 
figure 3 should be interchanged. 


New Table II 


TT 

p gm/cm 3 

/ 

8 

h 

R 2 

9mm ^ 9 ^ 9mn 
Condition 
Number 

^9^ 9mn I 
Condition 1 
Number || 

1 0.05 

0.0012 

3.0 

1 

0 

0.9936 

1.97 x 10 Y 

1.34 x 10 b 1 

0.10 

0.0012 

3.0 

1 

0 

0.9966 

2.92 x 10 s 

3.23 x 10 b | 

0.20 

0.0012 

3.0 

1 

0 

0.9933 

1.92 x 10 8 

2.12 x 10 b 1 

0.40 

0.0012 

3.0 

1 

0 

0.9997 

1.09 x 10 7 

1.73 x 10 b I 

0.60 

0.0012 

3.0 

1 

0 

0.9718 

6.54 x 10 8 

9.51 x 10 7 | 
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ABSTRACT 


A review of hot-wire anemometry was made to present examples of past work done in 
the field and to describe some of the recent and important developments in this extensive and 
ever expanding field. The review considered the flow regimes and flow fields in which 
measurements were made, including both mean flow and fluctuating measurements. 
Examples of hot-wire measurements made in the various flow regimes and flow fields are 
presented. Comments are made concerning the constant current and constant temperature 
anemometers generally in use and the recently developed constant voltage anemometer. 
Examples of hot-wire data obtained to substantiate theoretical results are presented. Some 
results are presented to compare hot-wire data with results obtained using other techniques. 
The review was limited to wires mounted normal to the flow in non-mixing gases. 


NOMENCLATURE 


°l ~ a 4 

constants in equation (47) 

a 

speed of sound 

A,B 

constants in equation (24) 

A, -A, 

constants in equation (18) 

b,-^ 

order of q in equation (47) 

B,-B } 

constants in equation (19) 

A'(Tj 

Lk, 

K 

overheat parameter, |(t^ogF w /^og/) 

B'{T 0 ) 

2Ljjtk,c f r w 

c . 

specific heat at constant pressure 

c. 

specific heat at constant volume 

C w 

specific heat of wire 

d 

wire diameter of mesh 


d( )fdt rate of change of quantity ( ) with respect to 
time 

d t diameter of cylinder 

d, diameter of jet 

</„ diameter of wire 

e' instantaneous voltage across the wire 
E mean voltage across the wire 

anemometer output voltage 
E' finite-circuit parameter, (l - s)/( 1 + 2 A^s) 
f frequency 

F dimensionless frequency 

F„ true one-dimensional spectral density 

F u measured one-dimensional spectral density 

turbulence reduction factor 


^ Formerly Senior Research Engineer, Analytical Services a Materials, Inc., Hampton, VA 33666 


£ 5 : - a 


Gr Grashof Number 

h coefficient of heat transfer 

h w height above wire shock generator to probe 

height above shock generator, immediate 

postshock value 
1 current 

k thermal conductivity of air evaluated at 

subscript temperature 
/r, wave number in the flow direction 
Kn Knudsen number 

L characteristic length 

m mean mass flow 

dlogp,/c\ogT, 

Mach number 
mesh size 

n exponent for mass flow in equations 

(16) & (17) 
n, ftogkJftogT, 

Nu Nusselt number evaluated at subscript 

temperature 

p mean static pressure 

p c mean total pressure 

P electrical power to the hot-wire 

Pr Prandtl number 

q sensitivity ratio, S m /S^ 

dynamic pressure 

Q forced convective heat transfer 

r sensitivity ratio, S m /S^ 

r d radial distance in cylindrical polar 

co-ordinate 

r p distance of virtual source of jet from origin 

r w radius of wire 

R resistance 

Re Reynolds number based on viscosity 

evaluated at subscript temperature and 
wire diameter 

mass flow - total temperature correlation 
coefficient, m’T p /mf v 

R, % velocity - total temperature correlation 

coefficient, u'T'Juf 0 

R^ density - total temperature correlation 

coefficient, ffTjpf 0 

R mp velocity - density correlation coefficient, 

u'ff/up 

R„ normalized auto-correlation function 

a sensitivity ratio, S fi /S^ 


S 


t 

T 

T f 

u,v,w 


“r 

X 

X. 
X „ 

y 

a 


P 

A 

S 

£ 

e 

% 


V 

9 

0 X 


A 

P 

r 

p 

T 

r w ' 

T ~aU 




sensitivity of hot-wire to the subscript 

variable 

time 

temperature 
(T.+ TJ/2 

velocity in x , y and z directions 

respectively 

frictional velocity 

distance in the flow direction 


virtual origin of the wake 
distance along the length of wire 
distance normal to the flow direction 




linear temperature - resistance coefficient 
of wire 
cb- \)M 2 


second degree temperature - resistance 
coefficient of wire 
boundary layer thickness 
displacement thickness for Blasius flow 
finite circuit factor, -(^og/./^ogR.^ 
transformed co-ordinate distance normal to 


body 

recovery temperature ratio, 7^/7^ 
temperature parameter, T„/T 0 
angle between plane sound wave and axis of 
probe 

mean free path 
absolute viscosity 
specific heat ratio, c p /c v 
density 
time lag 

temperature loading parameter, (r„ - T aA )/T Q 
temperature parameter, (T w - T^/T^ 
shear stress at the wall 
normalized fluctuation voltage ratio ,(%)/■$,. 


Subscript 

adw adiabatic wall condition 

adw,c adiabatic wall temperature, continuum flow 

condition 

adw,/ adiabatic wall temperature, free molecular 
flow condition 

B due to buoyancy effect 

C constant current anemometer 

e edge condition 
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eff effective velocity 

/ film condition 

h M and Re, are constant as Q varied 

0 total condition 

ref reference condition 

s electrical system untouched as M and Re, 

varied 

s sound source 

sc settling chamber 

1 evaluated at total temperature 

T constant temperate anemometer 

w wire condition 

ao free stream or static condition 

Superscript 

' instantaneous 

RMS 
mean 


INTRODUCTION 

Comte-Bellot 1 noted that the precise origin 
of hot-wire anemometiy cannot be accurately 
determined. One of the earlier studies of heat 
transfer from a heated wire was made by 
Boussinesq^ in 1905. The results obtained by 
Boussinesq was extended by King^ and he 
attempted to experimentally verify his theoretical 
results. These earlier investigations of hot-wire 
anemometiy considered only the mean heat transfer 
characteristics from heated wires. The first 
quantitative measurements of fluctuations in 
subsonic incompressible flows were made in 1929 
by Diyden and Kuethe 4 using constant current 
anemometiy where the frequency response of the 
wire was extended by the use of a compensating 
amplifier. In 1934 Ziegler^ developed a constant 
temperature anemometer for measuring 
fluctuations by using a feedback amplifier to 
maintain a constant wire temperature up to a given 
frequency. 

In the 1950's, Kovasznay®’^ extended hot- 
wire anemometiy to compressible flows where it 
was found experimentally that in supersonic flow 
the heated wire was sensitive only to mass flow and 
total temperature. Kovasznay developed a graphical 
technique to obtain these fluctuations, which is 
mostfy used in supersonic flow. In subsonic 


compressible flows the heat transfer from a wire is a 
function of velocity, density, total temperature, and 
wire temperature. Because of this complexity, these 
flow regimes were largely bypassed until the 1970's 
and 1980's when attempts were made to develop 
methods applicable^ for these flows. In recent 
years there were several new and promising 
developments in hot-wire anemometiy that can be 
attributed to advances in electronics, data 
acquisition/ reduction methods and new 
developments in basic anemometiy techniques. 

Previous reviews, survey reports, and 
conference proceedings on hot-wire anemometiy are 
included in references 1,9-20. Several books^l‘24 
have been published on hot-wire anemometiy and 
chapters^S'SO have been included in books where 
the general subject matter was related to 
anemometiy. 

This review considers the development of 
hot-wire anemometiy from the earliest 

consideration of heat transfer from heated wires to 
the present Although mean flow measurements 
are considered, the major portion of the review 
addresses the measurement of fluctuation 

quantities. Examples of some of the more important 
studies are addressed for wires mounted normal to 
the flow in non-mixing gases. The present review 
attempts to bring the development of hot-wire 
anemometiy up to date and note some of the 
important, recent developments in this extensive 
and ever expanding field. 

FLOW REGIMES AND FLOW FIELDS 

Based on the applicable heat transfer laws 
and suitable approximations, hot-wire anemometiy 
can be conveniently divided, into the following flow 
regimes: 

1. Subsonic incompressible flow 

2. Subsonic compressible, transonic, and 
low supersonic flows 

3. High supersonic and hypersonic flows 

Within each of these major flow regimes are the 
following sub- regimes: 

1. Continuum flow 

2. Slip flow 
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3. Free molecular flow 

In subsonic incompressible flow the heat 
transfer from a wire is a function of mass flow, total 
temperature and wire temperature. Since density 
variations are assumed to be zero, the mass flow 
variations reduce to velocity changes only. The non- 
dimensional heat transfer parameter, the Nusselt 
number, is usually assumed to be a function of 
Reynolds and Prandtl numbers and under most 
flow conditions the Prandtl number is constant. 
Evidence exist which indicate that Nu , is also a 

function of a temperature parameter®. In subsonic 
compressible, transonic and low supersonic flows 
the effects of compressibility influence the heat 
transfer from a wire. For these conditions the heat 
transfer from the wire is a f(u,p,T e ,T„) and 
Nu, = f(Re,,M,0). In high supersonic and 
hypersonic flows a strong shock occurs ahead of the 
wire and the heat transfer from the wire is 
influenced by subsonic flow downstream of the 
shock. Because of this, it was found experimentally 
that Nu, = f{Re„6) only, and the heat transfer from 

the wire is again a function of mass flow, total 
temperature, and wire temperature. 

In continuum flow the mean free path of the 
particles is veiy much less than the diameter of the 
wire and conventional heat transfer theories are 
applicable. Where the diameter of the wire 
approaches a few mean free paths between the 
particles, the flow does not behave as a continuum, 
but exhibits some effects of the finite spacing 
between the particles. These effects have been 
studied 31 * 3 ^ by assuming a finite velocity and a 
temperature jump at the surface of a body. This 
gas rarefaction regime was noted as slip flow. In 
free molecular flow the fluid is assumed to be 
composed of individual particles and the distance 
between the particles is sufficiently large that their 
impact with and reflection from a body is assumed 
to occur without interaction between the particles. 
Free molecular flow is theoretically studied 33 using 
the concepts of kinetic theory 3 **. 

Figure 1 presents a plot of Mach number 
vs. Reynolds number for lines of constant Knudsen 
number where «/. = 0.00015 inch and for flow 
conditions where 1.5 5 p a , psia % 150. Baldwin's 
results 33 indicated that the continuum flow regime 
existed for Kn < 0.001 and slip flow conditions 


existed for 0.001 <. Kn s. 2.0. Other references 
suggest that slip flow conditions were attained only 
for Kn > 0.01. Even using the larger value of Kn for 
the slip flow boundary, (i.e., Kn > 0.01) operating a 
0.00015 inch wire at low Mach numbers and at 
atmospheric conditions is near the slip flow 
boundary. If the total pressure is decreased or the 
wire diameter reduced, the value of Kn would be 
shifted farther into the slip flow regime. Free 
molecular flow conditions are generally assumed to 
exist for Kn > 2.0. Figure 1 can be used to delineate 
approximate values for M and Re, for the various 
sub- regimes. 

Various applications of hot-wire 
anemometry and the approximate level of velocity 
functuations are: 

Tvnes of Flows Approximate 39 

u ju 

0.05% 

0.20 - 2.00 % 

3.0 - 20.0 % 

2.0 - 5.0 % 
Over 20.0 % 


The two types of anemometers primarily 
used are the constant current anemometer (CCA) 
and the constant temperature anemometer (CTA). A 
constant voltage anemometer (CVA) is presently 
under development 57 . Even though these three 
anemometers are described as maintaining a given 
variable " consUmf , none of these strictly accomplish 
this. The degree of non-constancy for the CCA is 
determined by the finite impedance of its circuit 33 . 
The constancy of the mean wire temperature for a 
CTA at high frequencies is limited by the rate at 
which the feedback amplifier cam detect and 
respond to rapid fluctuations in the flow. The CVA 
maintains the voltage across the wire and leads 
constant rather than across the wire 57 . The non- 
constancy effects in the CCA and the CVA can be 


1. Freestream of wind tunnels 33 * 37 

2. Down stream of screens and 
grids 35 * 33 

3. Boundary layers 39 ' 43 

4. Wakes 43 * 44 

5. Jets 45 * 47 

6. Flow downstream of shocks 43 

7. Flight in Atmosphere 49 

8. Rotating Machinery 30 * 51 

9. Miscellaneous 53 ' 53 


TYPES OF ANEMOMETERS 
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accounted for by calibration of the CCA 88 and by 
knowing the lead resistance in the CVA. 

The heat balance for an electrically heated 
wire, neglecting conduction and radiation is: 

Heat Stored m Electrical Power In - Aerodynamic 


Heat Transfer Out 


— T w = P-Q 
dt " * 

(1) 

I 2 ^- iLdjfc-Tj) 

(2) 


If the heat storage term is properly compensated, 
then equation (2) becomes: 

/ 2 R,= xLdXT.-T^) (3) 

The measurement of fluctuations in a flow 
requires a sensor, in this case a wire, with a time 
response up to a sufficiently high frequency. The 
time constant of even small wires are limited and 
the amplitude response of these wires at higher 
frequencies decreases with frequency. Therefore, 
some type of compensation must be made for the 
wire output There are two methods for 
accomplishing this. Earlier approaches utilized a 
constant current anemometer with a compensating 
amplifier that had an increase in gain as the 
frequency increased 4 . An example of the roll off in 
the frequency of the wire, the gain of the amplifier 
and the resulting signal is shown schematically in 
figure 2. In principle, the output from the wire can 
be compensated to infinite frequencies. However, as 
the frequency increases, the noise output from the 
compensating amplifier will equal and ultimately 
exceed the wire output, which limits the gain that 
can be obtained. A schematic diagram of a CCA is 
presented in figure 3. 

The constant temperature anemometer uses 
a feedback amplifier to maintain the average wire 
temperature and wire resistance constant {i.e., 
dT,/dt = 0 in equation (2)}, within the capability of 
the amplifier. The practical upper frequency limit 
for a CTA is the frequency at which the feedback 
amplifier becomes unstable. A schematic diagram of 
a CTA is presented in figure 4. A third anemometer, 
presently under development 8 ^, is the constant 
voltage anemometer. This anemometer is based on 


the alterations of an operational amplifier circuit 
and does not have a bridge circuit A schematic 
diagram of a CVA is presented in figure 5. 

The upper frequency response of a CCA is 
generally accepted to be higher than that of a CTA. 
There is some evidence that the frequency response 
of the CVA might equal or exceed that of the CCA. 
The fluctuation diagram technique described by 
Kovasznay is usually used with a CCA to obtain 
data at supersonic speeds. This technique depends 
on the sensitivity of the wire being a function of wire 
temperature or overheat and the frequency 
response of the wire being assessable to 
compensation to almost zero overheat. This 
technique has limited application for a CTA, since at 
low overheats, the frequency response of the 
anemometer approaches the frequency response of 
the wire^ 8 . 

An example of the difference between the 
fluctuation diagrams obtained 89 using a CTA and a 
CCA is presented in figure 6. The intersection of the 
diagram with the vertical axis at = 0 

represents the total temperature fluctuation and the 
data show that the CTA cannot be used to measure 
these fluctuations. The reason for this is illustrated 
in figure 7a where the total temperature spectra at 
low overheats for the two anemometers are 
presented. In these cases the spectrum obtained 
with the CTA was attenuated at a frequency that 
was about two orders of magnitude less than for the 
CCA. At high overheats the two mass flow spectra 
were more nearly equivalent (figure 7b). However, 
in reference 60 the output of a laser was modulated 
and used to heat a wire to check the frequency 
response of a CTA. It was shown that the frequency 
response was essentially unchanged down to an 
overheat of 0.07. 

The CTA can be used to make 
measurements in supersonic flows ty using two 
wires. For these flows the CTA is operated with two 
wires having different but high overheats, digitizing 
the voltages and using two equations to obtain m ' , 
T' c and m'T e as a function of time 8 *. Then 
statistical techniques are used to obtain quantities 
of interest. In general, the CTA is more suitable for 
measuring higher levels of fluctuations than a 
CCA^ 8 . It remains to be determined how the CVA 
will compare with the CCA and CTA. At present it 
appears that the CVA has a higher signal to noise 
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ratio than either CCA or CTA. Additional 
advantages and disadvantages of the CCA versus 
the CTA are described in references 1, 29, 57, 62 
and 63. 

At low speeds a linearizer is often used to 
convert the non-linear relationship between wire 
voltage and velocity to a linear relationship. There 
are two types of linearizers in use; the logarithmic 
and the polynomial. A linearizer makes it possible 
to directly relate the measured voltage to the 
velocity. However, the linearization process does 
not result in better measured quantities^. 

LIMITATIONS OP HOT-WIRE ANEMOMETRY 

Most of the data obtained using hot-wire 
anemometry is limited to small perturbations. 
There are cases, however, where this linearization of 
the anemometry equation is not accurate and non- 
linear effects can influence both the mean^? and 
fluctuating^® voltages. Since high level fluctuations 
can influence the mean voltage measured across 
the heated wire, it is important to calibrate probes 
in flows with low levels of fluctuations. 

Because of the mass associated with the 
wire supports, there can be a significant amount of 
heat loss from the wire due to conduction to the 
relatively cold supports. This heat loss results in a 
!; panwise temperature distribution along the wire 
that, in turn, causes a variation of heat transfer 
from the wire^.21,64 along its length. In order to 
compare the heat transfer results from one wire or 
probe with another, the heat transfer rates must be 
corrected for these losses. However, computation of 
fluctuation quantities requires that the uncorrected 
values of the heat transfer rates be used. An 
example of the temperature distribution along a 
wire and its mean temperature^ * is shown in figure 
8. The finite length of the wire and its attendant 
temperature and heat transfer distribution 
influences the level of the spectra (especially at 
higher frequencies), correlations, and phase 
relationships between sensors^®’®®. 

The spacial resolution of a wire is limited by 
the length of the wire and the size of the smallest 
scale of fluctuations in the flow. If the length of the 
wire is larger than the smallest scale, the resultant 
magnitude of the spectra will be attenuated at the 


higher frequencies. The length of the wire with 
respect to the size of turbulence can have an effect 
on the measurements of fluctuation intensity, space 
and time correlations, and the turbulence scales 
and micro scales®®'®®. Additional spacial 
resolution problems encountered near walls was 
discussed in references 69 and 70. Proximity to 
walls of wind tunnels or to surfaces of models can 
introduce errors in measurements due to increased 
heat transfer from the wire due to conduction to the 
relatively cold walls^l»25 example of the effect 
of wire length on normalized spectra is presented in 
figure 9. The spacial resolution of multi-wire probes 
is further limited by the distance between the wires. 
The hot-wire probe intrusion into the flow can cause 
severe disturbance in certain flows. Examples are 
flows with large gradient such as boundary layers 
and vortices. Because of the above, hot-wire 
anemometry has limited resolution in space, time, 
and amplitude^. 

A severe problem is encountered in 
hypersonic flows when the gas is air. At higher 
Mach numbers the total temperature must be high 
enough to prevent liquefaction of air in the test 
section. There is a maximum recommended 
operating temperature for each wire material. 
These two facts places severe limitations on the 
maximum overheat at which wires can be operated. 
For example, the maximum recommended operating 
temperature for Platinum- 10% Rhodium wire is 
1842°R. For a M = 8 wind tunnel, the total 
temperature required can be as high as 1360°R. 
Using a recovery temperature ratio of 0.96, the 
maximum values for r„ is 0.394 and 9 W = 1.354. 
If gas rarefaction effects are experienced and rj is 
greater than one, then the problem is even more 
severe. For i) = 1.1 the maximum value for r w 
under the above conditions is 0.254. The above 
values for r. are based on the average temperature 
for the wire. For small L/rf. wires the limitation on 
r„ would be smaller due to higher temperatures at 
the mid-portion of the wire. The total temperature 
at low pressures where rj could be larger need not 
be as high as those at higher pressures, however, 
the constraint of constant total temperature during 
the calibration process limits the amount that T 0 
can be reduced. (Also see ref. 71-76). 
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PROBE PRE-CALIBRATION PROCEDURE 

Once a probe is constructed, the following 
procedure should ensure accurate and reliable 
measurements. First, the probe should be operated 
at the maximum 4 . and T w that will be used during 
the proposed test This is done to pre-stress and 
pre-heat the wire to ensure that no additional strain 
will be imposed on the wire during the test that 
could alter its resistance. For supersonic and high 
4 «, subsonic flows, the wires should also be checked 
for strain gaging, that is, stresses generated in the 
wire due to its vibration. Note, for testing in flows 
having high values of 4 ,, the wires should have 
slack to reduce the stress in the wires and to help 
eliminate strain gaging. If strain gaging is 
significant the wire should be replaced. During this 
pre-testing many wires might fail due to faulty wires 
or manufacturing techniques, but it is better that 
the wires fail in pre-testing rather than during an 
actual test 

A temperature- resistance relationship for 
wires is usually requires to compute the heat 
transfer rate from the heated wires. It is generally 
recommended that the following equation, which is 
a second degree equation in AT , be used: 

= 1 + ~ rj + Afc - (4) 

'v 

After the wires have been pre-stressed and pre- 
heated, they should be placed in an "oven" and the 
wires calibrated to determine the values for a, and 
/?, . Once this calibration has been completed, the 
probes can be placed in a facility for mean flow 
calibration over the appropriate ranges of velocity, 
density, total temperature and wire temperature. 

STATISTICAL QUANTITIES 

Data obtained using hot-wire anemometry 
are typically reduced to statistical quantities. Over 
the past few years the analysis of random data has 
been developed to a veiy high degree 77 ' 7 ®. This 
plus the rapid developments in electronics (i.e., the 
A/D converters and high speed computers), have 
made it possible to obtain almost any statistical 
quantity of interest within the error constraints of 


the heated wire. Much of this is due to the fact that 
the digital processing of data can be used to obtain 
many quantities that are difficult or impossible to 
obtain using analog data reduction techniques. 

Many types of single point and multi-point 
statistical quantities can be obtained using hot-wire 
anemometry® 0 '® 5 . j t j s routine to measure mean 
flow and RMS values, histograms and the higher 
order moments of skewness and kurtosis, auto 
correlation, and one dimensional spectra. 
Measurements of multi-point statistical quantities 
include cross correlations, two- point histograms 
and higher order two-point moments, cross spectra, 
and coherence functions. Attempts were made to 
measure higher moments up to eighth order® 1 . 

These measurements can be used in 
various ways to evaluate many characteristics of the 
flow such as scales, decay rate, energy content 
etc^ 5 . The coherence function is a useful statistical 
quantity that can be used to evaluate various 
properties of a flow® 4 . It can often be used to 
determine the predominant sound propagating 
angle and to determine the dominant mode present 
in a fluctuating flow field 85 *® 6 . 

A few examples of statistical quantities that 
were measured using hot-wire anemometry are 
presented in figures 10-13. Integral and micro time 
and length scales of a flow can be determined from 
autocorrelation functions such as the one presented 
in figure 10. The higher moments of skewness and 
kurtosis (figure lla-b) can be used to determine if 
the fluctuations are Gaussian. For a Gaussian 
distribution the value of the skewness parameter is 
zero and for the kurtosis the value is 3. Figure 1 1 
shows that both of these moments indicate that the 
mass flow and total temperature fluctuations are 
Gaussian over most of the thickness of the 
boundaiy layer. The value of third order auto- 
correlation function, such as the one shown in 
figure 12 , can be used to support turbulent flow 
theories. An example of space-time correlations 
measured in a turbulent boundary 5 ®*® 7 is 
presented in figure 13. The peak of these 
correlations at t * 0 indicate the presence of 
convection. The calculation of the convection 
velocity, obtained tty dividing the separation 
distance by the time at which the individual curves 
peaks, indicates that there was no significant 
variation of the convective velocity over the spacings 
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used. An example of normalized spectra measured 
downstream of a grid® 8 is presented in figure 14 
and show the increased attenuation of high 
frequency disturbances with increased distance 
downstream from the grid. (Also see ref. 89-92). 


MEAN FLOW MEASUREMENTS 

SUBSONIC INCOMPRESSIBLE - CONTINUUM FLOW 


GENERAL HEAT TRANSFER RELATIONSHIPS 

The heat transfer from a wire under the 
limits of the present report (i.e., the wires mounted 
normal to flow in non-mixing gases) is 2 ^ : 

Q= f(u,p t p,c f ,T w ,T^) (5) 

if the fluid properties of p , c f , and k are based on 
T c , then the above equation becomes: 

Q=f(u,p,T„A) (6) 

Since T^ = rjT a and ij= f(Kn,M)= f(u,p,T 0 ). For 
incompressible continuum flows equation (6) 
reduces to: 


Theoretical Considerations 

The functional relationship between the 
power to the wire or the heat transfer from the wire 
and the mean flow variables are required to 
determine the so called " static f calibration of the 
wire from which the sensitivities to the various flow 
variables can be obtained in order to calculate the 
fluctuations. Because of this the mean flow results 
and probe mean flow calibration procedure are 
considered together. 

The first attempt to obtain a theoretical 
solution for the heat transfer from a heated wire 
mounted normal to the flow was carried out by 
Boussinesq 2 . The equation that he obtained is: 

Q = L(2jxk,c,pur w )(T w - T^) (1 1) 


q - f(tt\ f t ) (7) Equation (1 1) can be expressed in terms of non- 

dimensional quantities as follows: 

Unless noted, the total temperature will be used 
throughout this report to evaluate p , , c r , and k, 

where as p will be based on T m . 



For a wire with a given L/d. the Nusselt 
number can be expressed 28 in terms of other 
dimensionless parameters as: 


Nu, = f\ 


T - T 
Re,,Pr,Gr,— ” 


t. 


( 8 ) 


King re-analyzed the problem of heat transfer from 
a heated wire and obtained the following 
relationship: 

Q = 4 *. +2j’*,c p pur w ){T w - O (13) 

or in terms of non-dimensional quantities: 


and can be written as follows to show the effects of 
compressibility: 

Nu, = /{he, ,Pr,Gr,M . (9) 

For relatively constant temperatures, Pr = constant 
and if Gr<Re *, buoyancy effects will be small and 
Gr can be neglected. These approximations lead to: 

Nu, = f(Re„M.,r w ) (10) 



From equation (11) and (13) it can be seen 
that the only difference between Boussinesq's and 
King’s results is the inclusion of the additional term 
k, in King's result that attempts to account for the 
effects of natural convection. At * high " values of 
Reynolds number the two results are essentially 
equal. 

Using equation (3), equation (13) can be 
expressed as: 
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(15) 


P = [AiT.)+BiT.UZ\[T w - T^] 


where A’{T 0 ) and B’{T 0 ) are based on King's results. 

However, the quantities A' and B' are usually 
determined for a given wire by direct mean flow 
calibration. Often the exponent for the mass flow 
term is determined from a curve fit to the data. The 
values for the exponent can ranged from 0.45 to 
0.50. 


For a CTA, equation (15) can be generalized 


to: 


where A r (Tj 


KMtJ 


= AriTj + BATjmf 

KBiL) 


«. Kr 

For a CCA: 


and B t {T„) = 


<*xKf 


(16) 


In this equation George et. al.,^8 noted that A, - A, 
are functions of T a . They proposed the following 
equation for the calibration of wires that is 
independent of T c for a limited range: 


Re, = B, + BjMr* + B 3 Nu +B t Nu* + B s Nu 2 (19) 
where p is evaluated at T 0 and k evaluated at T, . 

Example* of Data 

A summary of heat transfer data from 
cylinders in terms of Nu, vs. Re, taken in the 
subsonic continuum flow regime was presented in 
reference 21. The results from these experiments 
are compared with the theoretical results of 
Boussinesq and King in figure 16. This figure 
shows that there is a relatively good agreement 
between the measured results and King's theory 
over a wide range of Reynolds numbers. There is a 
substantial difference between Boussinesq's theory 
and King's theory and the measured results for 
Reynolds numbers less than about 100. 


uJ-o = 4c(r * ) + Bc(L){mT (17) 

where A c (tJ= A’(T.)/a,R^ and B C (T,)= B'(Tj/a t R«. 

Therefore, if wires operated with a CTA or a CCA are 
calibrated over a range of m and T 0 , equations (16) 
and (17) indicate that the calibration curves will be 
straight lines if the left hand side of the equation is 

plotted as a function of (mf . In general, the slopes, 
B, and intercepts. A, will be functions of the total 
temperature. If the total temperature is constant, A 
and B will be constants and if the density is 
constant, the mass flow term will reduce to velocity. 
An example of voltage versus velocity for a wire 
operated with a CTA is presented in figure 15 for 
various values of total temperature. (Also see ref. 
93-97). 

Hot-wires have also been calibrated in the 
form of u = /(is) rather than the more conventional 
form of E = /(«) . The constant T m and p version of 

King's law for a CTA is E 2 = A + Byfu and when 
expressed as u= f(E) gives: 

u = A, - AjE 2 + AyE 4 (18) 


A large amount of heat transfer data was 
also presented by McAdam" in terms Nu, vs. Re f 
and he recommended the following equation: 

Nu, = 0.32 + 0.43^ f- 52 (20) 

In comparison, King’s equation with Pr = 0.70 is: 

Nu, = 0.3183 + 0.6676(/fe,) a50 (21) 

For Reynolds number equal to zero the two 
equations give essentially the same value for the 
Nusselt number. At higher values of Reynolds 
number King's equation is about 40 percent higher 
than the values of Nusselt number presented by 
McAdam. A film temperature is often used as the 
temperature at which k , p and p are evaluated 
when correlating Nu versus Re. However, p is 
sometimes evaluated at the free stream static 
temperatures. The use of the film temperature for 
evaluating fluid properties has been questioned in 
reference 100. However, in this reference the 
density in the Reynolds number was evaluated at 
7), where as, in reference 99 this apparently was 

not the case. 
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Bradshaw^ notes that there is a difference 
of opinion throughout the hot-wire anemometiy 
community about the usefulness of a universal 
correlation based on variables evaluated at a film 
temperature. These correlations provide a useful 
guide for plotting results and comparing mean flow 
results obtained by different investigators. 
However, if good accuracy is to be obtained for the 
fluctuations, individual calibration of probes is 
required. 

Often attempts are made to measure mean 
velocities using hot-wire anemometiy. It can be 
shown using equation (13) that the voltage across a 
wire is a f(u,p,T t ,T w ). Therefore, to measure the 

mean velocity, the other variables must be held 
constant or a method must be used to correct the 
data for any variation in variables other than 
velocity 88 *^ *»*°*. Because of this, the hot-wire 
anemometer is not a very good mean flow 
measuring device, even if one utilizes some of the 
corrections that have been developed. However, for 
limited variations in the independent variables, 
corrected velocity measurements were reported in 
reference 101. 


Vary Low Valodtiaa 

At very low velocities the heated wire can 
cause a relatively significant vertical movement of a 
fluid due to buoyancy effects on the lower density 
fluid adjacent to the wire. This results in a change 
in the effective velocity around the wire. Efforts 
were made to calibrate and use hot-wire 
anemometiy at very low velocities 102-107 w here 
natural convection effects were present The 
influences of natural convection are parameterized 
by the Grashof number. Experimental evidence*® 8 
indicated that if Gr < Re 3 then free convection effects 
were negligible. An example of the effect of low 
velocities on the Nusselt number is presented in 
figure 17. 


King also provided an equation suitable for 
low speed flows: 


2nLk,(T„-T a J 

log (tyO 


( 22 ) 


where b- k,e'~ r /c f pu and y = Euler's Constant = 
0.57721. In terms of non-dimensional quantities, 
equation (22) becomes: 


Nu,= 


log[2e'-V(Pr*e,)] 


(23) 


Equation (22) is valid for ud w < 0.0187 where U is in 
cm/ sec and d. is in cm. Equation (13) is valid for 
ud . > 0.0187. 


For velocities as low as 1.0 cm/ sec, Haw 
and Foss* 0 ^ attempted to correlate their data using 
King’s equation in the form: 

E 2 =A + Bu n (24) 


A deviation of their data from a fitted curve was 
observed at u » 30 cm/sec. The diameter of the wire 
used in their experiment* 0 ^ was not noted. 
However, if one assumes a value of 0.00015 inch or 
0.00020 inch, the limits for the application of 
equations (13) and (22) indicate velocities of 49 
cm/sec or 37 cm/sec, which are not too different 
from 30 cm/ sec. The use of equation (22) would not 
improve the correlation presented in reference 102 
since it can be shown that as u -* 0 in equation (22) 
£„ -+ 0. The data of reference 102 indicates that at 
i# — 0 the intercept of the curve is greater than the 
value indicated by the intercept in equation (24). 
Correlations obtained using the results of a theory 
based on Oseen* 08 flow would not improve the 
correlation since this approach gives results that 
are similar to those obtained using equation (22). 
For a heated wire tested in horizontal wind tunnels, 
Utf cannot reach zero since the effective velocity is: 

u 2 =u 2 + u 2 (25) 

and for a heated wire u t * 0 . 


SUBSONIC SLIP FLOW AND TRANSONIC FLOW 


Theoretical Considerations 

These two flow regimes will be treated 
together since the experimental results are similar. 
Kovasznay 8 extended hot-wire anemometiy results 
to compressible flows and showed that there was a 
significant difference between the heat transfer in 
compressible and incompressible flows. Several 
experimenters obtained heat transfer 
measurements at low speeds and found an 
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apparent compressibility or Mach number 
effect®®’®^* at Mach numbers as low as 0.1. 
Spangenberg 1 conducted extensive tests over a 
wide range of variables and determined that the 
apparent compressible flow effects at Mach number 
as low as 0.05 was really due to gas rarefaction 
(e.g., slip flow). 

In this flow regime the heat transfer from 
the heated wire is generally given as: 

Q=P= nLk,{T„- jiT,)Nu, (26) 

In transonic flow and subsonic slip flows 
the Nusselt number is no longer only a function of 
Reynolds number and Kings' law is no longer 
applicable. The most common functional 
relationship for the Nusselt number in these flow 
regimes is®®: 

Nu, = f(M,Re„0) (27) 

since it was found that Nu, is also a function of a 
temperature parameter. Another functional 
relationship that was used to analyze gas 
rarefaction effects is®®: 

Nu, = r.) (28) 

In subsonic compressible flows the recovery 
temperature of the wire can change and functional 
relationships for tj are: 

»7 = f(M,Re,) or t,= f{M,Kn) (29) 



Dependent 

Variable 

Independent 

Variable 

Reference 

© 

Nu, 

R e ,> 0 

Morkovin® 8 


Nu, 

Kn, A/„, r. 

Baldwin 3 ® 

o 

Nu, 

Rc t , M m t 


0 

Nu, 

Kn, M„, 9 


© 

*1 

Re,, M. 

Morkovin® 8 

0 

n 

Kn, M m 

Vrebalovich 111 


Table I. Functional Relationships for Nu, and rj . 

Morkovin chose Nu, = f(M„,Re„0) and 
t)= f(M„,Re,) for the development of his equations. 
In order to emphasize the gas rarefaction effects. 


Baldwin chose Nu, = and 

tj = Independent variables that might be 

used to relate Nu, and rj to the dependent variables 
are presented in Table 1. Although Morkovin and 
Baldwin chose the variables of © and 0, one could 
just as well have chosen the variables noted in ® or 
0. It will be shown later that the variables in 0 
might be the most efficient group to use. 

Examples of Data 

Baldwin®® and Spangenberg**® 
investigated the heat transfer from wires over a wide 
range of Re,, M, and T w in the slip flow and 
transonic flow regimes. Their results, presented in 
figures 18a-b, shows that Nu, = f(Re,,M ) for Mach 

numbers ranging from 0.05 to 0.90 and Reynolds 
numbers ranging from 1 to about 100. The effects 
of wire overheat on the values of Nu, were also 
determined by Baldwin and Spangenberg and 
examples of these effects are shown in figure 19. 
The values of Nu, can increase or decrease with 
increased overheat depending on the Mach number 
and Knudsen number. 

Results from theoretical calculations made 
for the effects of slip flow on heat transfer from 
wires were reported in reference 31. An example of 
these results is presented in figure 20. The levels of 
the calculated Nusselt number do not agree with 
measured results, however, the trends of the 
theoretical results agree with the experimental 
trends shown in figure 18. 

SUPERSONIC CONTINUUM FLOW 

General results for compressible flow shows 
that Nu, = f(M,Re,,6). However, it was 

experimentally® determined that Nu, * f(M ) 

for Mach numbers greater than about 1.4. Typical 
heat transfer data for supersonic flow is presented 
in figure 21 to illustrate the approximate invariance 
of Nu, with M. At higher Mach number and 
relatively low total pressures, there is a high 
probability that much of the data presented for 
supersonic and hypersonic flows are in the slip flow 
regime. 
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FREE MOLECULAR FLOW 


Standler, Goodwin and Creager 33 
computed the heat transfer from wires for free 
molecular flow and an example of their results along 
with measurements are presented in figure 22. A 
combination of continuum flow, slip flow and free 
molecular flow results^® are shown in figure 23. 
From this figure it can be seen that for continuum 
flow at large Reynolds number Nu, * Re * . For free 
molecular flow Nu, « Re, and slip flow results 
smoothly connect the two regimes. Therefore, for 
slip flows, Nu, varies with exponent of Reynolds 
number which range from ‘/a to 1. 

RECOVERY TEMPERATURE RATIO 


The recovery temperature ratio must be 
known to compute the heat transfer from heated 
wires. In general, the recovery temperature ratio is 
a function of Mach and Reynolds numbers or Mach 
and Knudsen numbers. However, for Mach number 
greater than about 1.4 the recovery temperature 
ratio is not a function of Mach number for 
continuum flow. A " universaT curve presented by 
Vrebalovich 111 (figure 24) correlated the 
temperature recovery ratio with Knudsen number 
for all Mach numbers. Using the results presented 
in figure 24, the temperature recovery ratio for 
continuum flow and free molecular flow, curves of 
f} vs. M and Kn can be calculated. An example of 
these calculations is presented in figure 25. 


FLUCTUATION MEASUREMENTS 


SUBSONIC INCOMPRESSIBLE FLOW 

Theoretical Considerations 
a. Constant Temperature Anemometer 

For a constant temperature anemometer, 
King’s 3 equation can be expressed as: 

Tfl 

~ = *[*, + V2**,c,/*d.Jr„ - rjT c ] (30) 

where and T w are constants. 


If one assumes that the changes in k , , c p 
and rj can be neglected, the change in E will be a 
function of pu and T 0 as given by the following 
equation for small perturbations: 



( 31 ) 


where 




1 JlijRe,Pr 

4 \ + Jli^Re,Pr 


and S ^ 


1 7 

2 r. 


(32) 


From the above equation it can be seen that 
S m — ► 0 as Re, -+ 0 and S m — ► !4 as Re, — ► « . For the 
temperature sensitivity, — ► -oo as r.->0 and 

S K 0 as r„-*oo. Equation (31) shows that 
E = where S m = c\ogE/e7logm and 

S r< = <?log Efdiog T a . 

Since equation (31) shows that E = f(m,T c ), 
the fluctuation of mass flow and total temperature 
can be measured using a CTA 113 ' 113 . This can 
best be done by using two wires operated at 
different, but high overheats, digitizing the data, 
and solving two equations for m' , T p and m'7; as 
functions of time. 


If the total temperature and the Mach 
number varies significantly, then k, and c f must be 
differentiated with respect to T c and rj differentiated 
with respect to Mach number. Under these 
conditions it would be more appropriate to use the 
equation obtained by Rose and McDaid 8 with the 
assumption that Nu, * f{M). 


Instead of using King's equation, consider 
equation (15) for measuring mass flow and total 
temperature fluctuation. For a CTA equation (15) 
becomes: 


dioeB* 




r</ 10 g« + 


alr.X")"] 

, s_. 

A\T,)+ B\T.\m) n A\T.)+ B\T.lm) n 


A\T.) ^ A Y‘ ] (n.fBiT. fl 8 ™ 

cto 8 7, CTOgJ, rj 


L/logr. (33) 
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<note= 


y*°g— 


b. Constant Current Anemometer 

For the CCA anemometer, Kings' equation 
becomes: 


BirJmT 






El = L[k, + pnk,c f pud . |T. - r,T.] 

(34) 


[£r.)+ i-(r.X»)"j [£r> eiTjmf 


kh %T. 


(M) 


Again assume that k,, c f and rj are constant, the 
change in £ is given by the following equation: 


£ 

E 



(35) 


where 


«. *,r.(l-«) & 'xjRe.Pr 

' ^ ” [0 - - tfj (36) 


If dlogl - 0 then 

^ 'i r n J Re - Pr ... and S. = - k '- T L 




— ■ •: ■ fcUlU = - = 

2[k,T,-e] l + j2n^Re,Pr ‘ [k,r„~ 0} 


(37) 


Equation (36) and (37) shows that S m = /(r w ). If 
Re, -* 0 then S m -* 0 , but if Re,-* so then 




k,r„ 


2[k,r,-0] 


. On the other hand, if r„ -* 0 then 


S m -* 0 and if r w -* oo then S, 

If r, -> 0 then S T -* 

» r. Q 


2[l + V^V^]' 

and if r. -* oo then S T -* 0 . 


Again, it is possible to measure both m , f 0 and R mX> 

using a CCA 113 and the fluctuation diagram 
developed by Kovasznay^. An example of 
fluctuation diagrams for two discrete frequencies 
measured with a CCA is presented in figure 26. 
Again if the total temperature and the Mach number 
varies significantly, then it would be more 
appropriate to use Morkovin's equation with the 
assumption that Nu, * /(m). 


Similarly, to measure mass flow and total 
temperature fluctuation, equation (15) for CCA 
becomes: 


Mass flow fluctuations measured in 
subsonic flows can be very misleading where there 
is a significant amount of far- field sound. The mean 
square value of mass flow fluctuation is: 



The magnitude of the mass flow fluctuation depends 
on v , p , and R tip where -1 £ R tp £ 1 . As an example, 

assume R^,= 1 , indicating downstream moving 
sound, and u/u = p/p. Under these assumption the 
mass flow fluctuation equals twice the velocity or 
density fluctuation. However, if R up = -1, indicating 
upstream moving sound, and u/u = p/p; the mass 
flow fluctuation are zero. 

Examples of Data 

Most of the measurements made using hot- 
wire anemometry were and still are being made in 
the subsonic, incompressible, continuum flow 
regime. An extensive amount of data was 
accumulated over the years in various flow fields^ 3 . 

a. Freestream 

Some of the first fluctuation measurements 
made using hot-wire anemometers were obtained in 
the freestream of wind tunnels to help evaluate the 
effects of turbulence on the transition of laminar 
boundary layers * 26. The purpose of this effort was 
sm attempt to extend wind tunnel transition data to 
flight conditions in order that the on-set of 
transition might be predicted on full scale aircraft. 
Measurements in the freestream are also required 
to study the effect of freestream disturbances on 
laminar boundary layer receptivity. An example of 
measurements made in the freestream is presented 
in figure 27 for the Low Turbulence Pressure Tunnel 
located at the NASA Langley Research Center 3 ^. 
The filled symbols represent data taken in the 
facility during 1940^' and the curves represent 
measurements made in 1980. The agreement 
between the two sets of data is veiy good when it is 
noted that the data taken in 1940 was obtained at 
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p c = 4 atmospheres and the low datum point at 
Re, = 5 x 10 1 is for M » 0.02 . 

Fluctuation measurements were also made 
in various location within wind tunnel circuits, 
predominantly in the settling chamber. 
Anemometiy was used to evaluate the efficiency of 
contractions in reducing vorticity levels in the test 
section of wind tunnels* 18 . An example of the 
results obtained through a contraction is presented 
in figure 28. The absolute value of the velocity 
fluctuation in the direction of the flow was reduced 
through the contraction but the relative values were 
greatly reduced depending on the area ratio of the 
contraction. For example in figure 28 the velocity 
fluctuation downstream of the contraction are 
ratioed to the mean flow in the large section of the 
contraction where the local velocity is low. If these 
downstream velocity fluctuations were ratioed to the 
local mean velocity, these normalized fluctuations 
would be substantially smaller, i.e., uju H = 2.6 vs. 
«/« = 0.16. 

b. Grids 

It was found that screens or grids can 
effectively reduce vorticity fluctuations. Because of 
this attenuation, screens have been extensively 
investigated 25 ' 38,119 using hot-wire anemometiy 
to optimize their characteristics for use in wind 
tunnels to reduce the vorticity levels in the test 
section. An example of these measurements is 
presented in figure 29 where the turbulent 
reduction factor is given as a function of the A p/q. 
across the screens 12 ®. Therefore, the use of 
screens in the settling chamber along with a 
contraction of adequate area ratio, can substantially 
reduce velocity fluctuations in the test section due 
to voriticity 121 . (Also see ref. 122-130). 

c. Boundary Lavers 

Hot-wire measurements were made in 
turbulent boundary layers 42 to measure the 
Reynolds stresses and other fluctuation quantities 
to furnish data for the development of turbulent 
boundary layer theories. An example of 
measurements made in the boundary layer on a flat 
plate 131 is presented in figure 30a-b. Figure 30a 
shows the significant variation of the velocity 
fluctuation across the boundary layer while figure 
30b shows an example of the local stream wise 
velocity fluctuation ratioed to the local velocity. 
Forming the ratio in this latter manner indicates 


that the velocity fluctuations can exceed 40 percent, 
a value that is too large for an accurate assumption 
of small perturbation. 

Extensive measurements were made of 
turbulent flows in pipes 4 ®' 41 to compare 
theoretical and measured results. From these 
measurements many statistical quantities were 
obtained including Reynolds stresses; triple and 
quadruple correlations; energy spectra; rates of 
turbulent energy production, dissipation, and 
diffusion; and turbulent energy balance. An 
example of the streamwise velocity fluctuations 
across a pipe is presented in figure 31a-b. 

Hot-wire anemometiy has been extensively 
used to investigate the characteristic of various 
boundary layer flow manipulators such as Large 
Eddy Break-up devices (LEBUS) 132 , Riblets 133 and 
roughness elements 134 . Laminar boundary layer 
transition due to T-S waves 135 , cross flow 135 and 
Gortler vortices 13 ^ was extensively studied using 
the hot-wire techniques. Also the effects of heat 
addition 138 , sound 139 and vorticity 14 ® on 
boundary layer characteristics have been 
investigated. 

The hot-wire anemometer with a single wire 
cannot determine the direction of flow. However, a 
technique using a multi-wire, " ladder probe? was 
developed 141 to study the separated boundary 
layer where a significant amount of reverse flow 
occurred. This technique was used to determine 
the location of the zero average velocity in a 
subsonic turbulent boundary layer. (Also see ref. 
142-144). 

d. Laminar Boundary Laver Receptivity 

One of the major impediments to a through 
understanding of laminar boundary layer transition 
is the ability to predict the process ty which 
freestream disturbances are assimilated into the 
boundary layer. These free stream disturbances 
can be either vorticity, entropy, sound or a 
combination of these fluctuations. 

The effect of freestream fluctuations on the 
stability of the laminar boundary can be 
investigated by making measurements in the 
freestream and in the boundary layer to evaluate 
the receptivity of the boundary layer to 
fluctuations 145 in the freestream. An example of 
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fluctuations measured in a subsonic boundary layer 
on a flat plate for various frequency bands is 
presented in figure 32. Kendall * 45 presented three 
types of measurements made in a laminar 
boundary layer due to velocity fluctuations from the 
free stream. The first type is illustrated by the x's 
and consisted of broadband velocity fluctuation 
where the peak level occurs towards the inner part 
of the boundary layer. This type of measurement is 
noted as the KlebanofTs mode and is represented by 
the solid line. The results obtained when the data 
were filtered at the Tollmein-Schlichting (T-S) 
frequency are represented by circles. Although the 
frequencies were identical to those of T-S waves 
they were not T-S waves since the convection 
velocity was equal to the free stream value. The 
maximum level of these fluctuations occurred at the 
outer part of the boundary layer. The third type of 
fluctuation is represented by the dotted line. These 
were true T-S waves which occurred in packets and 
had a convective velocity of 0.35 to 0.4 of the free 
stream velocity. These peak fluctuation levels 
occurred near the wall. (Also see ref. 146-156). 

e. Jets 

Hot-wire measurements were obtained in 
jets 45 ' 47 to m easure the Reynolds stresses 
associated with free shear layers and to help 
evaluate the RMS levels and frequencies associated 
with jet noise. An example of the velocity and 
temperature fluctuations measured ^ 57 across a 
heated jet is presented in figure 33. The two types 
of fluctuations were normalized by the maximum 
and the local mean values, respectively. (Also see 
ref. 158-160). 

f. Wakes 

Various statistical quantities were 
measured downstream of a heated cylinder by 
Townsend 43 to obtain experimental results to help 
improve turbulent theories applicable to this type of 
flow. Some of the quantities obtained included 
turbulent intensities, sheer stress, velocity- 
temperature correlation, triple velocity correlation, 
diffusion rate and energy dissipation. 
Measurements were made from 500 to 950 
diameter downstream of the cylinder where 
dynamical similarity was assumed to exist An 
example of the mean curve fitted to the u- 
component of the velocity fluctuation is presented in 
figure 34 and shows good similarity. Uberoi and 
Freymuth 44 made extensive spectral 


measurements downstream of a cylinder and then- 
data indicated that only the spectra of large-scale 
turbulence were dynamically similar. 

SUBSONIC SLIP FLOW AND TRAN8QNIC FLOW 

Theoretical Considerations 

In compressible flows the heat transfer from 
a wire is usually described by the following 
equation: 

Q= «Lk,{T„- vTjNu, (40) 

Differentiating the above equation for the case 
where Q= P gives: 

diogP- —d logr. = rflogA'u, - — d\o$rj- -^-dlogT„ * dlogk, (41) 

T m T w T w 

The terms on the right hand side of the 
above equation depend only on the functional forms 
assumed for Nu, and tj (Table 1) and the chosen 
independent variables (Table II). Ultimately these 
terms depend on the variation of Nu, and 17 with the 
flow variables along with the aerodynamic and 
thermodynamic properties of the flow. The final 
form for the left hand side/ of the equation depends 
on the type of anemometer used. 

It was shown in reference 161 that, for a 
wire mounted normal to the flow, E = f{u,p,T.,T„) . 
Morkovin 58 and Baldwin 35 related Nu, and tj to 
the non-dimensional variables noted in Table I. 
However, recent results presented by Barre 
etal.,^ 5 ^ suggested advantages from using the 
following variables: E = and 

E = f{p,m,T 0 ,T w ). In order to obtain the equation for 
the CCA, they transformed Morkovin’s equations 
into their variables. This transformation is not 
necessary, since once the variables are chosen, the 
equations can be derived directly using a method 
similar to the one described by Anders 

The equations obtained using the variables 
of ® from table II gives the same results as those of 
O, Morkovin 58 , namely E - f(m,T 0 ) under the 

condition where Nu, * /( M m ) and = S p . Using the 
variables in ©, Barre etal., 15 ^ applied the 
resultant equations to measurements made in a 
turbulent boundary layer under the assumption 
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( 42 ) 


that p= 0 without assuming that S m = S p . They also 

p* u* ti 7' 

extended the equation to the case of supersonic flow _ = -S a — - S f! —+S i: — 

in the test section of wind tunnels by assuming that & u P T » 

all the fluctuations were sound, again not assuming 
S = S . where 

•t p 

Once the independent variables are chosen, 
it is not necessary to derive the equations using Nu, 
and t]. The " primitive? variables, u , p , T c etc., 
greatly simplifies the manipulation of the calibration 
data and can be used to correlate E as a function of 
» , p , T 0 etc., without evaluating Nu, and t ; . This 
technique might have advantages in the calibration 
of wires and the ease of operation of calibration 
facilities. 


J _ E'K r ( &o%Nu, ) flog A/m, 'j ( 1 <91ogl? <?logi ; i 

flo %u a <3ogAS J [a flog AS flog Re, j 


S _ (flogE _ E'Aj 


(flog Nu, (flog r] 
Oogp r„ (_ " cflog Re, cflogRe 


r„ 


-] 


(44) 


flog E ^ B*Aj 

r/ v 




flog Nu, 1 flogM/, 
flog/?#, 2a flog A/ 


-4 


1 < *° gly tm. a ° gn 
2a flog A/ ^flog/fr, 


:) 


( 45 ) 


The possible sets of variables based on the 
above discussion is presented in Table II. 



Dependent 

Variable 

Independent 

Variable 

Reference 

o 

E, Nu,, T] 

« ,P.T. 

Baldwin 35 , 

Morkovin 5 ® 

• 

E, Nu,, T] 

m , A/„, T a 

Barr, Quine and 
Dussauge 162 i 

B 

E, Nu„ rj 

P~. « . T c 

Barr, Quine and 
Dussauge 152 

1 

E 

« . P . T, 

Rose and 
McDaid® 
Stainback and 
Johnson® 5 

e 

E 

m . M m , T 0 


© 

E 

P„> » , T 0 



Table EL Various Independent Variables to 
Derive the Hot-Wire Anemometry 
Equations 

Therefore, there are many forms for the hot-wire 
equations depending on the variables chosen and 
the anemometer used. One should choose the 
variables that are most convenient for the flow 
situation under investigation. 


For a large range of Reynolds numbers and 
Mach numbers, 5. and S p in equation (42) are 

unequal 35 . Following Kovasznay's technique for 
supersonic flow, dividing equation (42) by , 

squaring and forming the mean gives: 




+ 2 ql R v ^-2qR ul ^-2sR. 
up *■ uT, 


'&■ ( 46 ) 

^pT. 


This is a general equation for a wire mounted 
normal to the flow in compressible flows where 
5. * S p . This is a single equation with six 

unknowns. In principal, this equation can be 
solved by operating a single wire at six overheats 
and solving six equations to obtain the three 
fluctuating quantities and their correlations. In the 
past, it was generally stated that the calibration of 
the wire cannot be made sufficiently accurate or the 
velocity and density sensitivities cannot be made 
sufficiently different to obtain a suitable solution 
using this technique. Demetriades 1 53 noted that 
the coefficient in equation (46) must occur to at 
least the fifth degree. This constraint, however, 
appears to be too restrictive. For example, assume 
that 4 is a function of q as follows: 

4=o,+ a i q t ' + a-#* 1 + arf* 1 (47) 


a. Constant Current Anemometer 

Using the heat transfer equation (40) and 
the functional relationship from equation (27) and 
(29), the change in voltage across a wire can be 
related to the changes in u , p , and T p . An example 
of a set of equations obtained for a constant current 
anemometer was given by Morkovin 5 ® as: 


It ca be shown that b , , b 2 and b } can have any value 
provided that the substitution of the relationship for 
4 into equation (46) results in an equation having 
at least six terms. An analysis of data obtained at 
transonic Mach number by S pangen berg indicates 
that 4 can be non-linearly related to q and suggest 
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HIGH SUPERSONIC AND HYPERSONIC FLOW 


b. 

Hors tm an and Rose 170 made 

measurements at transonic speeds where, for their 
flow condition, it was found that S, * S p . For this 

condition the transonic hot-wire problem 
degenerated to the supersonic flow problem where 
only m/m , fjT 0 and R kJ . could be measured. From 
their measurement of m/m , the velocity and density 
fluctuations were computed by assuming that fjT c 
and p/p were zero. An example of these results is 
presented in figure 40. In this figure Hors tm an and 
Rose's hot-wire results, represented by the circles, 
are compared with the thin film results obtained tty 
Mikulla 170 . 

c. Flight in Atmosphere 

Any attempts to extrapolate the effect of 
wind tunnel disturbances on laminar boundary 
layer transition to flight conditions requires some 
knowledge of the disturbance levels in the 
atmosphere. Much of the fluctuation data obtained 
in the atmosphere was measured using sonic 
anemometers on towers 179 . There was a limited 
amount of data obtained in the atmosphere using 
hot-wire anemometiy on flight vehicles 49 ' 180 . 
Otten et al., 49 expanded the methods devised by 
Rose and McDaid by using a two wire probe. One 
wire was operated by a CCA at a low over heat to 
measure f 0 . The other wire was operated with a 
CTA that was sensitive to m and f 0 . The results 
from these two wires were used to measure m and 

f e in the atmosphere. An example of spectral data 
obtained in the atmosphere is presented in figure 
41 and reveals the expected -5/3 slope, for m and 

t 


d. Subsonic Slip Flow 

For this regime Nu, = /(M„,Re,,r.) and 

S„ * S fi . These results are identical to those in the 
transonic flow regime and attempts have been made 
to apply the three wire technique developed for 
transonic flows to subsonic slip flows. For tests in 
subsonic slip flows the three wires were of different 
diameters in addition to being operated at different 
overheats. Some very preliminary data obtained 
using this technique in the Langley LTPT tunnel is 
presented in figure 42a-b where comparison with 
results obtained using King's equation are made. 


Theoretical Consideration 
a. Constant Current Anemometer 

In the 1950's and 1960's hot-wire 
anemometiy was extended into the high supersonic 
and hypersonic flow regime^’ 7 * 181,182 p or high 
supersonic flows it was found experimentally that 
S m = S p and equation (42) becomes: 



(49) 


Dividing equation (49) by the total temperature 
sensitivity, squaring, and then taking the mean 
results in the following equation: 



This equation was derived by Kovasznay 8 and used 
to generate fluctuation diagrams for supersonic 
flows. This equation was also used in references 
167 and 168 for subsonic compressible flows. The 
general form of equation (50) is a hyperbola where 
the intercept on the ^ -axis represents the total 
temperature fluctuation and the asymptotes 
represent the mass flow fluctuation 7 . 


Kovasznay demonstrated that the basic 
linear perturbation in compressible flows consists of 
vorticity, entropy and sound. He termed these basic 
fluctuations as " modes" . If the fluctuation diagram 
is assumed to consist of a single mode the diagrams 
were termed "mode diagrams''. An example of a 
general fluctuation diagram and the various mode 
diagrams for supersonic flow are presented in 
figures 43 and 44. 


b. Constant Temperature Anemometer 

For this case equation (42) becomes: 



(51) 


This is a single equation in two unknowns and a two 
wire probe can be used to obtain m ’ , TJ and m’T/ 

similar to the compressible subsonic flow 
0^61,183,184. 
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tomgl w of Data 
a. F r e es tr ea m 

In order to evaluate the relative 'goodness!' 
of supersonic wind tunnels and to relate the levels 
of disturbances in the test section to laminar 
boundary layer transition on models, a large 
amount of hot-wire measurements were made in the 
test sections of supersonic and hypersonic wind 
tunnels. In 1961 Laufer 188 presented 
measurements made in the test section of the Jet 
Propulsion Laboratory 18 x 20 inch supersonic 
wind tunnel over a Mach number range from 1.6 to 
5.0 using CCA. An example of the fluctuation 
diagrams obtained by Laufer is presented in figure 
45. Prom these diagrams the mass flow and total 
temperature fluctuations were obtained. Examples 
of the mass flow fluctuations are presented in figure 
46a. There was a significant increase of m/m with 
Mach number ranging from 0.07% at M = 1.6 to 
about 1.0 to 1.35% at M = S.O, depending on 
Reynolds number. All of the fluctuation diagrams 
were straight lines and Laufer demonstrated that 
these results indicated that the fluctuations were 
predominantly pressure fluctuations due to sound. 
Examples of the calculated pressure fluctuations 
are presented in figure 46b. Laufer concluded that 
the pressure fluctuations originated at the turbulent 
boundaiy on the wall of the tunnel and because of 
the finite value of the temperature fluctuations the 
sound source had a finite velocity. An example of 
the sound source velocities is presented in figure 
47. 

A large amount of hot-wire data was taken 
in the freestream of various facilities to measure 
disturbance levels in efforts to develop quiet 
supersonic wind tunnels. A review of this effort was 
reported in reference 186. 

Measurements in the freestream of the 
Langley Research Center Mach 20 High Reynolds 
number Helium Tunnel were performed by Wagner 
and Weinstein 18 *. All of their fluctuation diagrams 
were straight lines similar to the results obtained in 
supersonic flows. Examples of their measured 
mass flow and total temperature fluctuations are 
presented in figure 48. The mass flow fluctuations 
were substantially higher than the values measured 
tty Laufer at M = S.O . Pressure fluctuation 
measurements presented in figure 49 indicate that 
at low total pressures the boundaiy layer on the 
nozzle wall was probably transitional at the acoustic 


origin of the sound source. Relative sound source 
velocities are presented in figure 50. The source 
velocities for the Mach 20 tunnel at the higher 
pressures are significantly higher than those 
measured by Laufer at Mach numbers up to 5. 
(Also see ref. 187). 

b. Boundary Laver 

Measurements were made in supersonic 
and hypersonic turbulent boundaiy layers to extend 
the range of Reynolds stress measurements needed 
in the development of turbulent boundary layer 
theories. Barre eL al., 18 ^ conducted hot-wire tests 
in a supersonic boundaiy layer where transonic 
effects were accounted for by using a 
transformation of equation (42-45) from u , p , T c to 
p , m , and T a . Using the assumption that 
E = f{p ,m,T a ) and p/ p = 0, reduced their equation to 
E = f{m,T t ). Under these condition the fluctuation 
diagram developed by Kovasznay was used to 
obtain m/m, f 0 /T 0 and R m ^ without the assumption 
that S, = S p . 

Examples of their results are presented in 
figures 51 and 52. Figure 51 shows that the 
quantity ^/w' 2 /r. is greatly underestimated if the 
assumption is made that 5. = S p when the velocity 
in the boundaiy is transonic and S p . Figure 52 
show the variation of with the local Mach 

number through the boundary layer. The expected 
value for is -0.85 and the data obtained for 
* S p agrees well with this value. However, data 
evaluated where 5, was assumed to be equal to S p 
were substantially higher at the lower transonic 
Mach numbers. 

Fluctuations in a hypersonic boundaiy 
layer were made Ity Laderman and Demitrades 18 ^ 
and reported in reference 188. An example of the 
mass flow and total temperature fluctuations 
measured across the boundaiy layer is presented in 
figure 53. The velocity, density, static temperature, 
and pressure fluctuations were calculated using the 
mass flow and total temperature fluctuations and 
various assumptions. An example of these 
measurements is presented in figures 54 and 55. 

Additional measurements were made in 
hypersonic boundaiy layers by Laderman 18 ^ and 
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on an ogive cylinder by Owen and Horstman® 0 at 
M * 7.0 . The measurements made by Owen and 
Hors tm an included not only data for m, f e and i?.j- 
but integral scales and microscales, probability 
density distributions, skewness, kurtosis and 
intermittancy distribution across the boundary 
layer. A summary paper by Owen * 90 presents 
additional data which included space-time 
correlation, convective velocities, disturbance 
inclination angle, and turbulence life time 
distributions. (Also see ref. 87, 191-196). 

c. Wake 

An example of fluctuation diagrams 
measured in the wake behind a 15° half angle 
wedge at M = 1S.S is presented in figure 56. From 
these results Wagner and Weinstein concluded that 
the predominant fluctuation in this wake was 
entropy since the fluctuation diagrams were straight 
lines that intersected the r -axis at approximately 
-a. 

d. Downstream of Shock 

In supersonic and hypersonic flows the 
disturbances measured in the freestream of the test 
section are not necessarily the disturbances that 
can affect the transition of the laminar boundary 
layer on a body. The passage of sound waves 
through shocks will result in a combination of 
vorticity, entropy, and sound downstream of the 
shock 107 . Because of this the fluctuation diagram 
will no longer be a straight line but a general 
hyperbola^®. An example of this result is presented 
in figure 57. (Also see ref. 198-201). 

CONFIRMATION OF THEORETICAL RESULTS 

Hot-wire anemometiy was used extensively 
to validate or confirm theoretical results. Some 
examples of these efforts are presented below. 

Theoretical studies of the stability of 
laminar boundary layers to small disturbances were 
initially performed by Tollmein 202 ’ 20 ® and 
Schlichting 20 ^. fh ese calculations indicated that 
disturbances of a given frequency could decrease, 
remain constant or be amplified depending on the 
frequency chosen and the Reynolds number. The 
first experimental verification of the theoiy was 
made by Schubauer and Skramstad 1 ®®. An 

example of recent experimental and theoretical 


results for the determination of the neutral stability 
boundary 20 ® in a laminar boundary layer is 
presented in figure 58. 

Hot-wire measurements were made 
downstream of "grids" to evaluate the theory for the 
decay of turbulence. Tests conducted by Kistler 
and Vrebalovich 20 ® to evaluate the "linear" decay 
law is presented in figure 59 and confirm this law 
for large values of the Reynolds number. At lower 
Reynolds number 2 ® the exponent can be closer to 
1.20 - 1.25. Measurements of spectra for velocity 
fluctuations downstream of grids was also made 
and compared with theory. Two examples of these 
spectra are presented 20 ®’ 207 in figure 60a-b. The 
spectra in figure 60a had an insignificant amount of 
energy in the expected inertial sub-range indicated 
by a slope of -5/3. This result was attributed to the 
low Reynolds number of the flow. The spectra 
presented in figure 60b was measured in a high 
Reynolds number flow and revealed a significant 
amount of energy in the inertial sub-range. 
Theoretical calculations were made for the 
temperature spectra in a heated jet and 
downstream of a heated grid. An example of the 
theory and measurements 207 is presented in figure 
61. Attempts have been made to predict the 
influence of measured freestream fluctuations on 
laminar boundary transition. An example of this 
efforts 20 ® is presented in figures 62. 

A considerable amount of data was obtained 
by Stetson et al., 200 in hypersonic flow to study 
the stability of laminar boundary layer. An example 
of these results are presented in figure 63 and 
indicates the existence of first and second mode 
instabilities in the laminar boundary layer. These 
results were in agreement with those obtained by 
Kendall and Demetriades 2 1 0 and was in 
qualitative agreement with the theoretical results 
obtained by Mack*®®. 

OTHER APPLICATIONS 

Hot-wire anemometiy was used in shock 
tubes in an attempt to check the frequency 
response of probes and to trigger other 
events® 2 ’ 211 . 

Theoretical results indicated that there 
would be "temperature fronts or steps" in cryogenic 
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wind tunnels due to the injection of liquid nitrogen 
into the circuit Measurements were made using a 
hot-wire anemometer to determine the possible 
occurrences of these thermal steps 6 ^. 

A pulsed hot-wire was used to measure the 
velocities and flow angle in low speed flows^ 1 ^. A 
wire, which was operated by an anemometer, was 
place in the wake of a second wire which could be 
alternately heated. The time lag between heating 
the forward wire, this pulse being detected by the 
second wire and the distance between the wires was 
used to compute the velocity of the flow. In 
reference 213 a somewhat different technique was 
described that used two CCA's and a CTA to 
measure the velocity and flow angle. 

Hot-wire anemometry was used to obtain 
the location of transition in a laminar boundaiy in 
addition to obtaining some information on the 
fluctuations in the laminar, transitional and 
turbulent boundaiy layers 2 14,2 16 . in some flows 
where the fluctuation levels are high, such as a jet, 
a moving hot-wire probe was found to improve the 
accuracy of the results 216 . This technique is 
usually noted as flying hot-wire anemometry. (Also 
see ref. 217-221). 

Hot-wire anemometry was used 222 to 
measure the focal points of the laser beams for 
Laser Transit anemometry (LTA) Using a 

traversing mechanism and a CTA, the distance 
between the two beams was determined by 
measuring the difference between the two maximum 
voltage outputs from the anemometer. It was noted 
that additional information could be obtained such 
as the mean value of the beam intensity intersected 
by the wire, laser beam power, beam separation, 
beam diameter, beam divergence, cross-sectional 
beam-intensity distribution and relative beam 
intensity. 


CONDITIONAL SAMPLING 

Organized motion or structures in a 
turbulent boundary layer has been extensively 
studied 26,22 ^' 226 using the concept of conditional 
sampling. The flow of a turbulent boundary layer 
over a concave surface was studied in reference 226 
to search for organized motion in the boundaiy 


layer. The possible existence of organized motion is 
illustrated in figure 64 from measurements made 
with two hot-wire probes located 0.15 apart The 
conditional sampling technique was used to 
determine the characteristic shape of the mass flow 
signal during the passage of an organized motion. 
An example of these results are presented in figure 
64. Figure 65 shows that the measured event at an 
upstream and a downstream station had the same 
general characteristic shape. 

COMPARISON OF HOT-WIRE MEASUREMENTS 
WITH OTHER TECHNIQUES 

In the past the hot-wire anemometer, with 
all its limitations, was the only instrument available 
that was capable of measuring fluctuations with 
adequate fidelity 12 . To some extent, this is no 
longer the case as other techniques such as LV, LIF, 
CARS, and Raman scattering are now available for 
measuring various mean flow and fluctuating 
quantities. In inviscid flow where the fluctuations 
can be low, the anemometer is presently the only 
reliable instrument available. Compared to other 
techniques the anemometer is still relatively simple 
to operate and relatively inexpensive. Because of its 
long history, the results obtained using anemometry 
is still often used as a standard for evaluating 
measurements obtained using other techniques. 
The extent to which the anemometer cam maun tain 
these advamtages depends on the continued 
development of the other techniques. 

Tests were conducted in turbulent 
boundary layers to compare hot-wire results with 
other techniques to vadidate hot-wire/ laser 
velocimeter (LV) 22 ^> 22 ® amd hot-wire/ laser- 

induced fluorescence (LIF) 22 ^ techniques. An 
example of velocity fluctuations obtained using hot- 
wire anemometer and a LV system is presented in 
figure 66 for streamwise velocity fluctuations. The 
agreement between the two sets of data is very 
good. Measured density and static temperature 
fluctuations measured with a hot-wire amd a LIF 
system in a turbulent supersonic boundauy layer is 
presented in figure 67. Agaun, except for two values 
for the density fluctuation obtained with the LIF 
system, the agreement between the two sets of data 
is very good. Various experiments were made using 
LIF and LIF/RAMAN techniques to measure f t j T e 
amd pjp where results were compared with hot-wire 
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measurements 22 ®. An example 2 ®® of these results 
is presented in figure 68. The only disagreement 
between the hot-wire and the other results was 
attributed to a shock that apparently did not cross 
the hot-wire probe. 


CONCLUDING REMARKS 

A review was made to illustrate the 
versatility of hot-wire anemometry in addition to 
noting some of its limitations. The review included 
examples of results obtained in the various flow 
regimes and various types of flow fields. Examples 
of data were presented for the subsonic 
incompressible flow regime that were used to 
evaluate the flow quality in the test section of wind 
tunnels, to obtain measurements in turbulent 
boundary layer, and to substantiate or validate 
various theories of turbulence. 

Recently attempts to extend hot-wire 
anemometry into the transonic and subsonic slip 
flow regimes were presented for cases where S„ = S p 
and S u *S p . Examples of data obtained at high 
supersonic and hypersonic Mach numbers were 
presented. These results revealed that the 
fluctuation diagrams measured in the test section 
were a straight line, indicating that the 
disturbances in the test sections were due to sound 
mode. 

Examples were also presented to illustrate 
measurements made to substantiate turbulent 
theories, to compare with other techniques, and to 
illustrate the concept of conditional sampling. 
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Figure 1. Mach number vs. Reynolds number for 
lines of Constant Total Pressure and 
Constant Knudsen Number; 
d, = 0.00015 inch ; T Q = 560° F . 



Figure 2. Frequency Compensation for CCA. 



Figure 4. Schematic Representation of Constant 
Temperature Hot-Wire System; (ref. 62). 



Figure 5. Schematic Representation of Constant 
Voltage Hot-Wire System; (ref. 57). 
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Figure 6. Modal Analysis that Compares the CTA 
and CCA Systems; yjS= 0.23; (ref. 59). 


Figure 3. Schematic Representation of Constant 
Current Hot-Wire System; (ref. 62). 
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(a) Low Overheat; 

Total Temperature Fluctuations 
Figure 7. Normalized Power Spectra that Compare 
the Constant Temperature Anemometer 
(CTA) and Constant Current Anemometer 
(CCA) Systems; y/S= 0.23; (ref. 59). 



(b) High Overheat; 

Mass-Flow Fluctuations 
Figure 7. Concluded. 



Figure 8. Typical Temperature Distribution Along a 
Convection Controlled Hot-Wire; (ref. 21). 



Figure 9. Hot-Wire Spectral Transfer Function; (ref. 143) 



Figure 10. Auto-Correlation Coefficient R^ir); (ref. 25). 
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(a) Third Moment; Skewness 
Figure 11. Distribution of the Mass-Flow and Total 


Temperature Fluctuations Across the 
Boundary Layer; (ref. 59). 



(b) Fourth Moment; Flatness 
Figure 11. Concluded. 




Figure 13. Examples of Filtered Space-Time Correlation 
Coefficient - 4 kHz; (ref. 87). 



Figure 14. Decay of Energy of the Spectral 
Components; (ref. 88). 


Figure 12. Composite Third-Order Auto- Correlation 
Coefficient R e (r)tr, (ref. 25). 
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Figure 15. Voltage vs. Velocity for a Wire Operated 
with CTA; (ref. 93). 



Figure 16. Summary of Heat Loss from Circular 
Cylinders Over a Wide Range of 
Reynolds Number in Continuum Flow; 
(ref. 21). 



Figure 17. Interaction of Free and Forced 
Convection; (ref. 104). 



Figure 18a. Nusselt Number Correlation for 
Cylinders in Subsonic Slip Flow; 
T c = 80°F; 7. = 584'F; (ref. 35). 



Figure 18b. Spangenberg’s Reported Heat Loss 
Measurements from Electrically 
Heated Wires in Air; (ref. 171). 
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Figure 19. Normalized Nusselt Number Variation 
with Overheat for 0.00015 inch 
Diameter Platinum Wire Containing 
10% Rhodium; (ref. 110). 



Re, 

Figure 20. Predicted Nusselt Number Correlation 
from Approximate Slip-Flow Theory; 
(ref. 35). 



Figure 21. Summary of Supersonic Heat Transfer 
from Transverse Cylinders in Rarefied 
Air Flow; (ref. 21). 



Re, 


Figure 22. Comparison of Heat-Transfer Data with 
Free-Molecule-Flow Theory Using 
a= 0.57; (ref. 33). 



Re, 


Figure 23. Correlation of Convective Heat Transfer 
from Transverse Cylinders; (ref. 10). 
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Figure 24. Normalized Recovery Temperature Ratio vs. Figure 27. Velocity Fluctuations Measured in the 
Free Stream Knudsen Number; (ref. 21). Test Section of LaRC Low Turbulence 

Pressure Tunnel; (ref. 37). 
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Figure 25. Recovery Temperature Ratio vs. Mach Figure 28 Effect of 16:1 contraction on Turbulence 
Number for Constant Values of Generated by 2-inch Square Mesh Grid; 

Knudsen Number. R . 3710; (ref. 1 18). 
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Figure 26. Fluctuation-Diagram of the Filtered 
Signals (Single- Wire); (ref. 97). 
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Figure 29. Measured Turbulence Reduction Factor 
for Various Screens and Combinations 
of Screens with Honeycomb Compared 
With Different Theories; (ref. 120). 



(a) Zero Pressure Gradient 
Figure 30. Relative Turbulence Intensities in a 
Boundary Layer Along a Smooth Wall; 
(ref. 131). 



(b) Constant-Stress Layer 
Figure 30. Concluded. 





(a) Across the pipe 

Figure 31. Distribution of velocity fluctuation in a 
pipe; (ref. 40). 
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(b) Near the wall 
Figure 31. Concluded. 



Figure 32. Types of Fluctuation Simultaneously 
Present in Boundaiy Layer; (ref. 145). 
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Figure 34. Relative Turbulence Intensity (M- 
component) in the Wake of a Circular 
Cylinder; Re^ = 1360; (ref. 43). 



Figure 35. General Fluctuation diagram for Subsonic 
Compressible Flow; (ref. 168). 


Figure 33. Distribution of Relative Intensity of Turbulent 
Velocity and Temperature Fluctuations in a 
Hot Round Free Jet; (ref. 157). 
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(a) Disturbances from Single Sound Source 
Figure 36. Fluctuation Diagram; M= 0.76; (ref. 167). 
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(b) Combination of Vorticity and 
Entropy Modes Downstream of a 
Cylinder 

Figure 36. Concluded. 
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Figure 37. Nusselt Number as a Function of p ; (ref. 110). 



(a) Velocity Fluctuations; T 0 = 303° K ; r„ = 0.672 . 
Figure 38. Spangenbeig; Data; (ref. 163). 
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(b) Density Fluctuations; T c = 303° K; 
r. = 0.672 . 

Figure 38. Continued. 
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(c) Total Temperature Fluctuations 
T, = 303° K\ p= 0.0008 gm/cm 5 . 

: 38. Concluded. 
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(b) Density Fluctuations 


Figure 39. Continued. 


(a) Velocity Fluctuations 
Disturbances Measured in the LaRC 8' 
Transonic Pressure Tunnel using a 
Three Element Hot-Wire Probe; (ref. 86). 


(c) Total Temperature Fluctuations 
Figure 39. Continued. 
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Figure 41. Spectra of Mass Flux and Total 
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(a) Measured Fluctuations 

Figure 42. Fluctuations Measured in the LaRC Low 
Turbulence Pressure Tunnel; 
p 0 = 15 psia; (ref. 172). 
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Figure 44. Generalized Mode Diagram for 
Supersonic Flow. 


(b) Mass Flow Fluctuations 
Figure 42. Concluded. 
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Figure 43. Generalized Fluctuation Diagram for 
Supersonic Flow. 
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Figure 45. Mode Diagram Measured in Test Section 
for M- 5.0; (ref. 185). 
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(a) Mass Flow Fluctuations 
Figure 46. Variation with Tunnel Mach Number; 
(ref. 185). 
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(b) Pressure Fluctuation 
Figure 46. Concluded. 




Figure 47. Variation of Source Velocity with Tunnel 
Mach Number; (ref. 185). 
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Figure 48. Root-Mean-Square Mass-Flow and Total- 
Temperature Fluctuations in the LaRC 
22" Hypersonic Helium Tunnel; (ref. 181). 



Figure 49. Pressure Fluctuation in the LaRC 22" 
Hypersonic Helium Tunne; (ref. 181). 



Figure 50. Relative Source Velocities in the LaRC 22" 
Hypersonic Helium Tunnel; (ref. 181). 
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Figure 51. Velocity Fluctuation Measurements in 
Boundary Layer; Influence of the 
Sensitivity Coefficient to Mach 
Number; Shaded Zone - Range of the 
Available Velocily Fluctuation Data in 
Supersonic Boundary Layers, 
Including Hot-Wire and Laser-Doppler 
Measurements; (ref. 162). 
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Figure 53. Variation of Wide-Band RMS Mass-Flux 
and Total-Temperature Fluctuations 
Across the Boundary Layer. Kistler's 
Results at M = 4.76 (Solid Line) are 
Included for Comparison; (ref. 188). 
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Figure 54. Variation of Wide-Band RMS Velocily and 
Static Temperature Fluctuations Across 


Figure 52. Variation of with Mach Number the Boundary Layer; p' * 0; (ref. 188). 


Through Boundary Layer; (ref. 177). 
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Figure 55. Variation of Wide-Band RMS Pressure 
and Density Fluctuations Across the 
Boundary Layer; p'* 0; (ref. 188). 



Figure 57. Fluctuation Mode Plots; Me 21; 

p, = 12.41 jc 10 s N/m 2 ; (ref. 48). 
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Figure 56. Mode Diagrams in Wake of 15° Half- 
Angle Wedge; (ref. 181). 


Figure 58. Neutral Disturbance in Studies of 
Laminar Boundary Layer Stability; 
(ref. 205). 
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Figure 62. Effect of Free-Stream Turbulence on 
Boundaiy Layer Transition; (ref. 208). 
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Figure 63. Amplification Rate vs. Frequency, 
Comparison with Others; (ref. 209). 
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Figure 64. Typical Fluctuating Mass Flow Data 
from Two Hot-Wires Indicating 
Similarity Between the Signals; (ref. 
226). 
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Figure 65. Ensemble Averaged Positive Events in 
the Weak Perturbation Case; (ref. 226). 
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Figure 66. Comparision of Streamwise Velocity 
Fluctuations Measured by Hot-Wires 
and Laser Velocimeter; (ref. 227). 



Figure 67. A Comparison of the Distribution of RMS 
Fluctuation Amplitudes in Static 
Temperature and Density within the 
Boundary-Layer obtained using a Hot- 
Wire at Single and Multiple Overheat 
Ratios with Direct Measurements 
Obtained using Laser-Induced 
Fluorescence (LIF). The LIF data have 
been Corrected for Instrument Noise; 
(ref. 229). 



Figure 68. Comparisons of Temperature and 
Density RMS Fluctuation Amplitudes 
Obtained with LIF Techniques or 
Implied by Hot-Wire Anemometer 
Data. Flow Conditions are the Same 
as in Figure 66. For Each Variable, 
the Fluctuation Amplitude is 
Normalized by the Local Time- 
Averaged Value; (ref. 230). 


52 


8tainback, P.C. and Nagabushsna, KLA. 


